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I 


A  three  part  report  is  presented  describing  the  mode  competition  in  a  gyroton 
oscillator  with  two  competing  electromagentic  modes.  Part  I  gives  the  theory 
of  this  mode  coupling  is  presented  using  a  quasi-1 inear  assumption  which  is 
valid  for  small  power  levels.  The  results  are  interpreted  in  terms  of  the 
phase  plane.  These  general  results  are  applied  to  a  specific  case  of  mode 
competition,  that  between  TEjqj  and  TEqh  modes  in  a  cavity  of  rectang^ular 
cross  section.  Part  II  presents  the  theory  of  coupling  which  amy  occur 
between  any  nutnber  of  cavity  modes  through  finite  conductivity  in  the  cavity 
walls,  holes  in  the  conducting  cavity  walls,  or  though  interaction  with  an 
electron  beam.  Part  III  details  the  experimental  work.  The  electrodynamic 
circuit  is  described  as  are  microwave  and  electron  beam  diagnostics.  Results 
of  initial  tests  are  also  given. 
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Competition  Between  Electromagnetic  Modes 
in  a 

Free-Electron  Maser 

Abstract 

A  three  part  report  is  presei\teci  describing  the  mode  competition  in  a  gyrotron  oscillator 
with  two  competing  electromagnetic  modes.  Part  I  gives  the  theory  of  this  mode 
coupling  is  presented  using  a  quasi-linear  assumption  which  is  valid  for  small  power 
levels.  The  results  are  interpreted  in  terms  of  the  phase  plane.  These  general  results  are 
applied  to  a  specific  case  of  mode  competition,  that  between  TE]0i  and  TEqh  modes  in  a 
cavity  of  rectangular  cross  section.  Part  II  presents  the  theory  of  coupling  which  may 
occur  between  any  number  of  cavity  modes  through  finite  conductivity  in  the  cavity 
walls,  holes  in  the  conducting  cavity  walls,  or  though  interaction  with  an  electron  beam. 
Part  III  details  the  experimental  work.  The  electrodynamic  circuit  is  described  as  are 
microwave  and  electron  beam  diagnostics.  Results  of  initial  tests  are  also  given. 

This  work  has  been  published  in  part  in 

A.H.  McCurdy  and  J.S.  Plewa,  IEEE  Trans.  Plasma  Sci.  20, 139  (1992), 
and  has  been  presented  at  the  following  conferences; 

A.H.  McCurdy  and  R.  Liou,  ICOPS 1994,  Santa  Fe,  New  Mexico. 
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PARTI 


Mode  Interaction  through  Amplitudes  and  Phases 
in  a  Two-Mode  Gyrotron  Oscillator 

Abstract  -  Part  I 

An  analysis  is  made  of  the  temporal  evolution  of  electromagnetic  modes  in  a  two-mode  gyrotron 
oscillator  characterized  by  phase  and  amplitude  interaction  through  the  terms  linear  in  the  oscillator 
power.  The  problem  is  solved  in  the  context  of  amplitudes  and  phases  which  vary  slowly  compared  to 
the  period  of  oscillation.  Specific  reference  is  made  to  competition  between  TEnq  modes  in  a  closed 

cavity  gyrotron.  Qualitative  features  which  are  found  include  phase  locking,  beat  frequencies,  periodic 
pulling  and  mode  excitation.  This  work  has  applicability  when  the  frequency  separation  between  the 
modes  is  on  the  order  of  the  frequency  bandwidth  of  each  mode,  or  the  modes  are  equally  spaced  in 
frequency.  Gyrotrons  of  this  type  include  those  with  low  quality  factor  modes  or  degenerate  modes. 
Phase  interaction  in  the  case  of  equally  spaced  cavity  modes  is  of  importance  in  analyzing  mode  locking 
phenomena. 
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I.  INTRODUCTION 


Operation  of  high  power  microwave  oscillators  at  high  frequencies  typically  involves  the 
competition  between  many  electromagnetic  modes.  In  vircators  [1]  and  gyrotrons  [2],  researchers  have 
attempted  to  optimize  the  devices  for  power  and  efficiency  in  spite  of  multi-mode  effects.  In  many  cases 
the  modes  are  spaced  far  apart  in  frequency  compared  to  the  modal  linewidth.  If  the  modes  are 
unequally  spaced  in  frequency  then  they  will  interact  primarily  through  their  amplitudes.  However,  in 
the  cases  where  1)  the  modal  quality  factors  are  such  that  modes  are  approximately  equally  spaced  in 
frequency:  2  (O2  -  0)3  -  (Oj  ~  0)2/Q  (Q  is  the  loaded  cavity  quality  factor)  or  2)  the  modal  lineshapes 
directly  overlap  (O2  -  ©1  <  0)2  /Q  then  the  relative  phase  between  the  modes  will  play  an  integral  role  in 
the  temporal  evolution.  Several  factors  tend  to  enhance  the  probability  that  one  of  these  cases  will  be 
realized  in  practice.  First,  high  frequency  oscillators  typically  have  closely  spaced  modes  simply  because 
the  mode  spacing  varies  as  X,”,  where  X  is  the  wavelength  and  n  is  a  positive  number.  Secondly,  it  is 
often  desired  to  keep  the  energy  density  stored  in  the  cavity  at  a  minimum.  For  a  high  output  power 
oscillator  this  implies  that  the  cavity  modal  quality  factors  will  be  low.  Thus  it  is  important  to  study  the 
case  where  phase  interaction  between  modes  is  possible.  A  closely  related  problem  is  that  of  the  phase 
locking  of  two  strongly  coupled  oscillators,  such  as  magnetrons  [3].  It  is  to  be  noted,  however,  that  in  the 
magnetron  locking  problem  there  is  an  intrinsic  time  delay  due  to  the  physical  separation  of  the 
oscillators.  Here  the  modes  coirununicate  instantaneously  since  they  are  not  separated  in  space. 

This  work  complements  previous  theoretical  work  on  gyrotron  mode  competition  in  that  1)  it  is 
in  the  spirit  of  the  analytical  rate  equation  approach  of  [4]  rather  than  the  largely  numerical  work  of  [2], 
15],  and  [6]  and  2)  it  emphasizes  the  low  quality  factor  regime  where  phase  interaction  occurs.  The 
complication  due  to  the  inclusion  of  phases  in  the  rate  equations  goes  so  far  as  to  provide  phase  coupling 
through  the  linear  terms.  Hence  there  is  no  regime  of  unimpeded  linear  modal  growth.  Of  particular 
concern  here  is  the  temporal  evolution  of  the  modes  and  the  nature  of  the  final  steady  state.  An  analysis 
of  the  effects  of  these  phase  dependent  linear  terms  reveals  phenomena  similar  to  those  seen  when  an 
external  signal  is  injected  into  a  single-mode  oscillator.  One  practical  consequence  of  this  work  is  that 
two-mode  oscillators  may  well  produce  single  frequency  output,  since  the  modes  are  generally  phase 
locked. 

This  work  examines  the  coupling  of  two  modes  through  amplitudes  and  phases  and  makes 
specific  application  to  the  closed  cavity  gyrotron  oscillator.  The  rate  equations  for  the  electromagnetic 
modes  are  given  in  Section  II  and  the  linear  cross-coupling  coefficients  are  calculated  for  the  gyrotron. 
The  rate  equations  are  analyzed  in  the  linear  regime  in  Section  III  in  order  to  quantitatively  describe  the 
phase  locking  phenomenon.  Typical  parameters  are  given  for  the  gyrotron.  In  Section  IV,  the  effect  of 
the  nonlinear  rate  equation  terms  is  outlined  with  emphasis  on  phase  locking. 
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II.  GYROTRON  RATE  EQUATIONS 


Here  we  focus  on  the  sinrplest  case  of  phase-amplitude  mode  interaction,  that  of  two  modes  with 
an  arbitrary  frequency  separation.  A  restriction  will  be  made  to  only  odd  powers  of  the  field  amplitude 
in  the  electronic  susceptibility.  Microwave  devices,  such  as  the  gyrotron  oscillator,  utilizing  a  phase 
mixed  electron  beam  demonstrate  this  property. 


A.  Temporal  Evolution  of  Electromagnetic  Modes 

A  set  of  equations  for  the  evolution  of  the  jf*'  electromagnetic  mode  can  be  written  [7] 


^  +  c2  kj2  aj  (t)  COj  X  Lj  j-  ^  ay  (t)  =  -  c2  ^  d^r  Jose  .  e* 


(1) 


where  the  complex  modal  amplitudes,  electron  beam  current  contributions  and  cavity  quality  factors  are 


e*  =  ijj(t)ei[vj(t)  +  wjt] 


and 


(l+i)ki  ^ 

kj-  Qojj' 


1 


(2) 


respectively.  Here  Aj(t)  and  yj(t)  are  the  slowly  varying  amplitude  and  phase  of  the  cavity  electric  field 
and  Jj{t)  is  the  slowly  varying  compx}nent  of  the  electron  beam  current.  The  beam  current  can  be  related 
to  the  cavity  fields  by  the  electron  equations  of  motion 

where  Otj^j  are  the  linear  growth  and  cross  coupling  coefficients  and  Pj  is  the  coefficient  of  self-saturation. 
Higher  order  saturation  and  phase  interaction  through  the  nonlinear  terms  is  neglected  in  this  analysis. 
Using  (2)  and  (3)  in  (1),  equations  for  the  modal  evolution  are  found 
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(4) 


where  ttj  is  the  linear  growth  rate  of  the  jth  mode. 

The  modal  rate  equations  are  obtained  from  (4)  by  using  the  method  of  slowly  varying  amplitude 
and  phase  [4],  [7],  [81.  These  equations  are 


2  0=  (Oi  po  +  coi  a’l  -  0)2  QL2  +0)1  (^2/pj^[(pi2  +  a’i2 )  COS  O  +  (ai2  -  ot’i2 )  sin  o]  + 

0)2  [(^21  -  “21 )  sin  O  -  (p2i  +  aji )  cos  o]  . 


Here  the  over  dot  indicates  differentiation  with  respect  to  time.  The  second  order  (in  power)  cross¬ 
coupling  has  been  neglected  as  well  as  frequency  pulling  or  pushing  due  to  the  second  order  interaction. 
In  (5)  Pj  is  the  power  in  the  f  ^  mode  (proportional  to  Aj^)  and  primes  used  on  the  complex  coefficients 
indicate  the  real  and  imaginary  parts:  x  =  Xj  +  i  xj".  The  actual  frequency  of  oscillation  of  mode  j  is  0)j 
and  O  =  \|r2  -  Vi  +  (CO2  -  C0|)  t.  The  cavity  effects  are  taken  into  account  through: 


Po  =  — 

0), 


0)i  0)2 


1^11  ^22 
iQo  Qo 


0), 


P21  = 


0)o 


0)2  Qo 


P2I  “ 


^12 

“iQo 


Oij  =  -I-  —  and 


^jk  = 


kj(l  +i)  1  1 

0). 

Here  0)°  is  the  cold  (beam-free)  cavity  resonant  frequency  and  and  indicate  the  quality  factors 
due  to  Ohmic  losses  and  external  losses,  resjjectively.  The  superscripts  indicate  that  these  losses  can  be  a 
source  of  mode  coupling  (these  factors  can  be  derived  using  the  formalism  of  Slater  [9]). 


7 


In  order  to  make  these  equations  as  nwthematically  tractable  as  possible,  we  assume  that  the 
cross-coupling  due  to  the  cavity  walls  and  external  load  is  small.  Thus  and  »  Qq®  and 
This  allows  the  coefficients  Pjj^  and  to  be  neglected  in  (5). 

Because  of  the  neglect  of  terms  which  are  higher  powers  of  Pj  in  these  rate  equations  it  is  clear 
that  this  representation  will  have  validity  only  in  the  small  signal  regime.  Thus  the  phenomena  described 
here  will  have  applicability  to  the  growth  stage  of  high  power  oscillators.  However,  in  a  nonlinear 
oscillator,  the  final  steady  state  is  intimately  dependent  on  the  characteristics  of  oscillation  growth.  Thus 
it  is  expected  that  many  of  the  features  found  here  will  persist  into  the  steady  state. 


B.  Coefficients  of  the  Rate  Equation  Expansion 

In  order  to  calculate  the  electron  beam  dependent  coefficients  of  the  rate  equations,  an  averaging 
is  done  over  the  electron  motion  through  the  oscillator  cavity.  The  beam  contribution  to  the  mode  is 
obtained  from  the  electronic  susceptibility  of  that  mode  (as  defined  in  [7])  which  is 


Z  =  - 


iJi 


co  e  A. 

J  O  j 


In  the  case  of  a  closed  cavity  gyrotron  with  a  circular  cross-section,  the  susceptibility  can  be  written  in 
terms  of  the  electron  momentum  as  [10): 


iR^ 


o  . 
CO-  e  A- 
J  o  J 


(6) 


where  and  qj  are  the  azimuthal,  radial,  and  axial  mode  indices  for  mode  j,  fq  is  the  axial  electric 

field  profile  (the  asterisk  indicates  complex  conjugate  and  z  is  the  axial  coordinate),  L  is  the  cavity  length, 
and  pj^  and  A  are  the  slowly  varying  components  of  the  electron  perpendicular  momentum  familiar  from 
gyrotron  theory  [11].  The  electron  mass  is  mg,  initial  electron  velocity  perpendicular  and  parallel  to  the 
axial  direction  are  v^q  and  Vzq  =  Vj^g  /c).  Iq  is  the  dc  beam  current,  Rq  is  the  guiding  center 
position,  and  the  perpendicular  wavenumber  is  kj^/  is  a  normalization  constant  which  is  given  in 
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Ill].  The  brackets  indicate  an  average  over  initial  conditions  of  the  electron  beam:  <j)Q  is  the  phase  angle 
on  the  electron  gyro-orbit,  Gq  is  the  guiding  center  phase  angle,  and  tg  is  the  entrance  time  of  the  electron 
into  the  cavity.  From  the  electron  equations  of  motion  the  momentum  can  be  expanded  in  powers  of  the 
electric  field  amplitude: 


-iAj  -iA:Z  -i((Ojt,-<>J 

PiC  '=Pi.e  '  e  ' 


1  +  X  S  anplzlx^’lx*) 


n=l p=0 


where  the  a^^p  are  z  dependent  expansion  coefficients  and 


'lyr 


iAr 


For  two  modes  the  averaged  momentum  expansion  can  be  written: 


£  b„j(z) 
s  =  1  n  =0 


2 

k  =  l 


,2n 


(7) 


where  the  b^j  are  expansion  coefficients.  Keeping  only  the  n  =  0  and  1  terms  in  (7),  (this  includes  the 
third  order  terms  in  the  momentum  expansion),  averaging  over  guiding  center  angle,  and  substituting 
into  (6),  the  following  expression  is  obtained  for  the  susceptibility 

Xj  =  aj  +  ajk  ^  e-i  |3j  a} 


where  the  definitions  for  the  (complex)  linear  growth,  linear  coupling  and  self-saturation  terms  are 
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“j  =  hjCn,.^j  |fq.{zfbo.(z)dz  ^  ajk  =  hjC^4  6m.,n^  I  f*(z)fq,{z)bo,(z)d2 
/laO  Jt<"0 


and 


Pi  =  hi  [  If, 

/s«0 


(8) 


.{zf  bij(z)dz 


Here  hj  =  i  Rj  /  (Eq  COj).  The  coefficients  ttj  and  (Jj  have  been  given  elsewhere  [4,  7].  The 
coefficients  CXjj^  are  found  by  using 


boj(z)  = 


fqe‘Aj^dz*-i 


fqe*Ai^'  dz" 


in  (8).  Note  that  there  is  no  phase  coupling  through  the  linear  terms  if  the  modes  have  different 
azimuthal  indices  mj.  Numerical  examples  of  Ojk,  calculated  for  TEni  and  TE|]|2  modes  in  a  circular 

cylindrical  gyrotron  cavity  are  shown  in  Figs.  1  and  2.  The  electron  beam  parameters  are  20  kV,  1  A, 
perpendicular-to-parallel  velocity  ratio  of  1.4,  and  Rq  =  .9  cm.  The  cavity  parameters  are  L  =  60  cm,  wall 
radius  of  1.79  cm  and  quality  factor  of  120.  Figure  1  shows  the  real  and  imaginary  parts  of  the  TE|n 

linear  growth  rate  and  linear  cross-coupling  (X|2-  Tho  coefficients  are  shown  as  a  function  of 

2 

normalized  detuning  of  the  wave  frequency  from  the  cyclotron  frequency,  Ac ;  = - r —  (cO;  -  Q).  Here 

'  “oPio 

Q  is  the  electron  cyclotron  frequency  and  (Oq  is  a  normalization  frequency  (tOg  /  27C=  5.0  GHz  in  the 
example  shown).  It  can  be  seen  in  Fig.  1  that  the  cross-coupling  coefficients  are  generally  an  order  of 
magnitude  smaller  than  the  linear  growth  rates.  That  the  cross-coupling  acts  to  oppose  the  linear  growth 
can  be  seen  from  Fig.  1  where  the  0,12  (both  real  and  imaginary  parts)  are  180°  out  of  phase  with  a|. 
The  TE||2  results  show  similar  features  with  respect  to  phase  and  relative  size.  A  comparison  of  Fig. 
2(a)  with  Fig.  1(a)  shows  that  0,21  =  ■  0t|2.  This  relationship  holds  between  any  modes  with  adjacent 
axial  mode  indices.  It  can  further  be  shown,  for  modes  spaced  more  than  one  index  apart: 

(Ojc  ajcj  =  COj  ajic  if  j  and  k  are  both  even  or  both  odd, 

cok  ttkj  =  -  COj  ajk  if  either  j  or  k  (but  not  both)  is  odd.  (9) 
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These  symmetry  relations  only  hold  when  Ac  j  =  Ac  For  most  values  of  the  magnetic  field,  the 

cyclotron  detuning  of  the  two  modes  is  not  the  same.  The  cross-coupling  coefficients  are  given  as  a 
function  of  magnetic  field  in  Fig.  3  for  the  TEm  and  TE|  ]2  modes.  It  can  be  seen  that  there  is  a  slight 

off^t  between  each  set  of  curves.  Hence,  when  solving  (5)  for  a  typical  set  of  gyrotron  parameters,  the 
coupling  coefficients  will  not  obey  (9). 

Figure  4  shows  the  coupling  coefficient  between  TEjn  and  the  TEj  J3  and  TEi^4  modes.  The 
result  here  is  similar  to  that  of  coupling  with  the  TE|i2  ”^ode,  except  that  the  strength  is  somewhat 
weaker  and  the  coupling  coefficient  changes  sign  more  rapidly  with  a  change  in  A.  The  degree  of 
coupling  of  the  TE|n  with  higher  order  modes  is  shown  in  Fig.  5.  The  coupling  strength  falls  off  asq"^, 
where  q  is  the  axial  mode  index.  The  effect  of  higher  order  modes  on  the  TEjj  j  is  expected  to  be  small 

since  the  frequency  spacing  between  modes  increases  with  q.  Hence  these  modes  make  a  small,  rapidly 
varying  contribution  to  the  power  of  the  TEm. 

Equations  (5)  are  first  studied  in  the  case  where  the  pj  =  0.  In  this  case  the  mode  coupling  occurs 

only  through  the  terms  which  are  first  order  in  the  p)ower.  It  will  be  seen  that  the  overwhelming 
tendency  is  for  the  modes  to  phase  lock  (that  is  reach  a  situation  where  dC>/dt  =  0).  Since  only  two 
modes  are  considered  here,  we  note  that  d<I>/dt  =  0  implies  that  the  two  modes  have  the  same  final 
frequency,  thus  they  are  also  frequency  locked  (or  synchronized).  This  phase  locking  occurs  whether  or 
not  the  powers  Pj  reach  a  saturated  level.  Application  is  made  to  the  TEj  jj  -  TEn2  competition  in  the 

closed  cavity  gyrotron. 


m.  Small  SIGNAL  Analysis 


One  interesting  feature  of  this  problem  is  the  fact  that  coupling  can  occur  between  the  modes  due 
to  terms  directly  proportional  to  the  mode  power.  These  terms  can  thus  cause  a  "linear"  mode  interaction 
of  the  same  size  as  the  linear  growth  rate.  Since  most  oscillations  start  from  a  small  noise  level,  thus 
spending  a  fair  aiiK)unt  of  time  in  the  linear  regime,  it  is  of  interest  to  examine  the  linear  mode  coupling 
predictions. 

A.  Solution  of  the  Linearized  Rate  Equations 

Equations  (5),  with  Pj  =  0  and  in  normalized  form,  become 
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where  t=rit,A=CO|  (pg  +  Otn’  -  a.22  ©2/ 0)i)/ (2  F^),  P  =  P/F  j  ,  F  =  F  2/F^,  and 
C|  =  (Xp  /F I  ^  C2  =  OCj  /F|  where  Fj  =  COj  (Oj"  -  Oj)  and  £0j  0.\2  =  0^  +  i  (Xj.  It  is  assumed  that  the 
cross-coupling  coefficients  have  the  symmetry  of  the  TE-jn  -  TE112  gyrotron  modes  given  in  (9).  The 
variation  in  CCji^  due  to  the  unequal  cyclotron  detunings  of  the  two  modes  is  neglected  (this  correction  can 
be  included  at  the  cost  of  complicating  the  following  analysis).  A  is  the  difference  between  the  free- 
running  beam  and  cavity  loaded  resonant  frequencies  of  the  two  modes  (normalized  to  the  growth  rate  of 
mode  1). 

It  is  to  be  noted  that  when  Pj  =  0,  (4)  represent  a  linear  system  of  equations  which  can  be  solved 

exactly.  For  two  modes,  the  solutions  are  of  the  form 
4 

aj(t)=  X  Pkje^‘ 

k=i  (11) 

where  the  coefficients  pi^j  are  found  from  the  initial  conditions  and  the  are  the  eigenvalues.  The 
eigenvalues  are  found  from  a  fourth  order  characteristic  equation  (with  complex  coefficients),  which  can 
be  derived  from  (4)  using  standard  techniques.  Instead  of  solving  the  characteristic  equation  exactly,  we 
focus  on  some  qualitative  features  of  the  solutions.  In  particular,  from  (11)  it  is  clear  that  the  ratio  of 
power  levels  in  the  two  modes  generally  approaches  p^i  / Pn2/  where  is  the  eigenvalue  with  the 
largest  real  part.  If  Re{A^}  <  0  then  the  power  in  both  modes  approaches  zero  (mode  suppression),  and 
if  Re{Xj|}  =  0  and  lm{Xj^}  ^  0  then  the  modal  power  levels  oscillate  in  time. 

Equation  (10c)  )delds  a  phase  locking  bandwidth  if  the  power  level  ratio  of  the  two  modes  is 
fixed.  This  stationary  phase  result  (d^/ dx  =  0)  is  possible  when  the  following  criterion  is  met 


In  the  case  where  the  cavity  and  load  cross-coupling  have  not  been  neglected,  the  locking  criterion  is: 
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where 


— 2A 
V 


<  1 


(pl2*‘^l2)-{^'/pJ‘{p2l  +  “21) 

b=P^J‘(a;2-a:2)-P>4fe-a2,) 

This  criterion  is  an  extended  version  of  that  obtained  in  single  mode  injection  locked  oscillators  [12].  This 
equation  differs  from  that  of  [12]  in  that  the  two  power  levels  can  be  of  comparable  size  and  beam 
coupling  effects  (as  well  as  coupling  through  the  external  load)  are  included  in  the  coefficients  a  and  b. 
Equations  (11)  and  (12)  indicate  that  the  two  modes  will  generally  phase  lock  together  through  the  linear 
interaction.  To  evaluate  the  locking  bandwidth,  one  requires  knowledge  of  the  steady  state  modal  power 
ratio.  This  information  must  come  from  a  solution  of  (10). 


To  examine  the  phase  locking  phenomenon  it  is  convenient  to  make  the  transfonnation  (Pj  ,  P2  ) 
(X/ y)  where  X  =  (P|  /  P2  and  y  =  (Pi  P2  Equations  (10)  result  in 


^  =  |-[l  -  r  +  (^  -  x)  Cl  sin  <I)  +  (^  +  x)  C2  cos  <t>]  (a) 

^  ~  2 ^  (x  ~  x)  C2  cos  O  +  (i  +  x)  Cl  sin  o]  (W 

=  A  -  ^  -  x)  C2  sin  <I>  +  ^  +  x)  Cl  cos  C>  (c) 


The  advantage  of  this  transformation  is  that  a  decoupling  is  achieved  between  y  and  the  other  two 
variables.  Thus  the  variable  y  has  no  effect  on  the  dynamics  of  x  and  O.  In  addition  the  instabilities, 
inherent  in  a  system  with  no  high  order  saturation,  are  confined  to  the  variable  y.  The  phase  locking 
predictions  of  this  system  of  equations  will  now  be  discussed. 


B.  Linear  Phase  Locking 
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dx  do 

The  two-mode  system  reaches  a  phase  locked  equilibrium  when  both  and - are  zero.  From 

dt  dx 

(13a)  and  (13c)  this  happens  for  an  angle  O® 


sin  O 


Cl  (r-  i)  +  2AC2 

xj{Ci2  +  C22) 


cosO  = 


C2(r-  l)-2ACi 
(i+  X«)(C,2  +  C22) 


where  is  a  real  (positive)  solution  to  the  real  quartic  equation 


(14) 


-  aop3  +  ai  p2-  aop+  1  =0  where p  =  (xO)^and 


(15) 


^  =  ,  ,n([Cl(r-l)t2AC2Ft  [C2(r.l)-2AC,P|  , 

(Ci^  +  C2W 


ai  =  - - 2 — 

The  stability  of  these  equilibria  can  be  found  by  a  perturbation  analysis  of  (13a)  and  (13c).  It  is  found  that 
the  growth  (or  decay)  rate  of  perturbations  is  given  by 


([C2(r-l)-2AC,p-[Ci(r-l)  +  2Acj)-2  . 


rp= 


(i-r) 


cosO 


Cl  sin  0°] 


(16) 


Stable  equilibria  (Fp  <  0)  correspond  to  phase  locked  solutions  of  (10).  Note  that  this  locking  takes  place 
in  the  absence  of  any  nonlinearity.  It  is  of  interest  to  determine  the  time  required  for  the  modes  to  achieve 
the  phase  locked  state.  An  estimate  of  this  lock-in  time  can  be  made  from  (16).  The  real  part  of  the 
eigenvalue  resulting  from  the  perturbation  analysis  gives  the  rate  of  approach  to  the  phase  locked 
equilibrium  point  (in  the  immediate  vicinity  of  that  point).  However  this  rate  is  found  to  be  quite 
accurate  even  at  points  far  from  equilibrium.  The  lock-in  time  X  lcx:k  defined  as  the  time  required  for 
the  oscillator  to  decrease  its  distance  from  the  equilibrium  point  in  the  X  -  O  plane  by  a  factor  of  c  =  2.718 
and,  according  to  (16)  is: 
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'tLOCK=— ^ 
-Tp 


(17) 


If  the  phase  plane  distance  is  measured  as 


then  a  comparison  can  be  made  between  the  predicted  rate  of  approach  to  equilibrium  and  that  found 
numerically  from  (10).  This  comparison  is  shown  in  Fig.  (6)  for  parameters  corresponding  to  gyrotron 
operation  at  a  magnetic  field  of  1.838  kG  (Fj  =  0.756,  F  =  1.8,  C|  =  0.0574,  C2  =  -0.0249,  A  = 
-0.978).  The  lock-in  time  in  this  case  is  3.29  ns,  which  is  -  16  radiation  periods.  It  can  be  seen  that  the 
value  of  X  lock  ^17)  provides  an  excellent  approximation  to  the  actual  lock-in  time.  The  actual 
approach  to  the  phase  locked  state  involves  a  ^ight  oscillation  which  is  not  found  from  the  approximate 
method.  Thus  the  time  required  to  reach  the  phase  locked  state  can  be  accurately  predicted  even  though 
the  modal  power  levels  are  far  from  stationary.  This  is  an  important  result  since  high  power  oscillators 
are  operated  in  the  short  pulse  regime  (the  power  changes  rapidly  in  time)  and  the  time  available  for  the 
locking  process  is  restricted. 

Phase  locking  properties  are  calculated  for  gyrotron  operation  when  the  TE|n  mode  is  linearly 

unstable.  This  unstable  region  (Fig.  7)  is  between  1.83  and  1.86  kG,  using  the  operating  parameters  of  Fig. 
1.  It  can  be  seen  that  the  TE112  niode  is  linearly  unstable  at  low  magnetic  field  and  is  not  normally 
excited  above  1.845  kG  (in  the  region  a|"  >  0).  The  stationary  amplitude  ratio  and  phase  locking  angles 
are  shown  in  Fig.  8  for  the  gyrotron  through  this  magnetic  field  range.  This  figure  shows  that  there  is  a 
discontinuity  in  the  stationary  equilibrium  (x®,  C)®)  at  1.841  kG.  At  this  value  of  the  magnetic  field  F  =  1 
and  a  periodic  solution  results  (infinite  lock-in  time).  The  values  of  X®  on  either  side  of  the  discontinuity 
are  approximately  reciprocals  of  one  another.  Above  1.841  kG  the  TE||2  mode  is  excited  (though  it  has  a 
negative  growth  rate)  by  energy  influx  from  the  TEui  mode.  For  lower  magnetic  fields,  the  TE112 
mode  dominates,  and  linearly  suppresses  the  TEm  mode  (the  amplitude  ratio  in  this  region  is  about  x® 
=  .03). 

It  is  of  interest  to  find  the  dependence  of  X  lock  o"  parameters  A,  F,  C|  and  C2  .  This 
information  is  given  in  Fig.  9.  Here  the  parameters  are  varied  about  the  gyrotron  coefficients  given  in 
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connection  with  Fig.  6.  In  Fig.  9,  the  four  curves  show  the  result  of  varying  one  parameter  at  a  time.  The 
detuning  is  varied  in  Fig.  9(a).  Unlike  the  injection  locked  oscillator,  the  minimum  locking  time  is  slightly 
offset  from  A  =  0,  and  the  modes  do  not  break  out  of  the  phase  locked  condition  at  large  detuning. 
Instead,  X  lcX!K  asymptotes  to  t ««  =  2/ (F-l).  This  is  due  to  the  dependence  of  the  steady-state 
amplitude  ratio  on  A.  Since  1  /(x°)^  +  (x°)^  continually  increases  with  A,  the  locking  limits  given  by 
(12)  are  never  met  (for  F  ^  1).  The  results  in  Figs.  9(a)  and  9(c)  thus  show  lock-in  times  relative  to  X  oo- 
These  results  will  not  be  accurate  when  A  becomes  very  large  with  F,  and  C2  fixed.  This  is  because, 
from  (12),  in  order  to  have  a  phase  locked  solution  with  large  A  either  x®  »1  or  x®  «1.  In  either  case, 
the  power  in  one  of  the  modes  is  very  large,  thus  violating  the  linear  assumption  made  in  (10).  The 
asynunetry  of  the  locking  results  in  Fig.  9(a)  is  due  to  the  presence  of  both  sin<[)  and  COsO  components 
in  (13c).  This  equation  does  not  have  even  parity  in  A,  and  (13a)  and  (13b)  do  not  have  even  parity  in  O. 

From  (14),  (15),  and  (16)  it  can  be  shown  that  X  lock  depends  on  F  in  the  form  (F-1),  thus,  the 
relevant  physical  quantity  involving  growth  rate  is  the  difference  in  growth  rates  of  the  two  modes. 
Figure  9(b)  shows  that  for  small  cross-coupling,  the  lock-in  time  depends  on  I  F  -  1 1 .  The  time  required 
for  phase  locking  becomes  very  large  at  F  =  1  since  X  00  becomes  infinite  at  that  point.  When  F  =  1,  the 

lock-in  time  becomes 


'Clock= 


(Cl^4-C2^)[l-(x0rt 

4ACiC2(x0f 


(18) 


It  is  found  that  when  A  Cj  C2  >  0  then  x®  <  1,  and  when  A  Cj  C2  <  0,  then  x°  >  1.  In  this  way, 
”^LOCK  always  positive.  However,  when  A  Cj  C2  =  0,  then  a  periodic  solution  may  be  obtained.  This 
solution,  depending  on  the  initial  conditions,  may  be  either  bounded  or  unbounded  in  O.  The  former 
corresponds  to  periodic  mode  pulling,  while  the  latter  is  mode  beating.  Another  possibility,  when  A  Q 
C2  =  0/  is  that  X®  =  1,  thus  yielding  a  rinite  ratio  in  (18)  and  a  phase  locked  solution. 

Generally,  the  more  disparate  the  growth  rates,  then  the  shorter  the  lock-in  time.  This  result  is 
interpreted  as  a  donunation  of  the  mode  interaction  by  the  more  rapidly  growing  mode,  hence  causing 
more  rapid  phase  locking. 

The  dependence  of  X  loCK  o*'  ^2  shown  in  Fig.  9(c).  The  figure  shows  that  an 

increase  in  the  size  of  the  cross-coupling  coefficient  I  Cj  I  causes  the  modes  to  phase  lock  more  rapidly. 
The  opposite  is  true  for  IC2  I .  Thus,  the  larger  I  Re{ai2}  I  and  the  smaller  I  Iin{aj2)  I  the  more 
rapidly  will  phase  locking  be  obtained. 


IV.  NONLINEAR  Saturation  (p  O) 
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The  inclusion  of  self-  saturation  effects  in  (13)  causes  a  considerable  analytic  complication.  When 
P  0  the  dependent  variable  y  appears  in  (13a)  thus  modifying  the  x  -  <I>  dynamics.  Thus  the  problem  is 
fully  three  dimensional  and  the  phase  space  is  more  difficult  to  visualize  than  in  the  linear  case.  A 
general  overview  is  given  of  the  numerical  results  and  application  is  made  to  the  TEm  -  TE|i2 
competition  in  the  gyrotron.  Only  the  case  Pj/  P2  <  0  is  considered  since  it  is  that  case  which  provides  a 
stable  steady  state  for  the  two  modes. 

It  is  found  numerically  that,  even  with  very  small  self-saturation  coefficients  P,  the  phase  locking 
between  modes  can  be  broken  by  an  increase  in  the  detuning  when  F  >  0.  This  is  so  because  p  prevents 
either  modal  power  from  becoming  extremely  large,  and  the  linear  growth  rates  and  cross-coupling 
prevent  the  powers  from  becoming  extremely  small.  Thus  the  allowable  range  of  x°  values  is  restricted. 
From  (12),  it  is  seen  that  the  phase  locking  frequency  band  depends  on  1  /  (x°)^  -t-  (x°)^.  Because  of 
the  restriction  on  x°,  the  denomirutor  of  (12)  cannot  increase  without  limit  to  balance  a  large  numerator 
(as  in  the  linear  case).  It  is  found  that  the  steady  state  values  x°  and  <I>°  depend  of  the  ratio  P1/P2  and 
are  thus  virtually  independent  of  P  in  the  case  P|  =  P2.  The  actual  values  of  the  steady  state  modal 
power  levels  are  sensitive  to  the  value  of  P,  however  this  information  is  confined  to  the  variable  y. 

When  P  ^  0  and  F  >  0  a  progression  of  three  oscillator  behaviors  is  seen  as  A  is  increased 
through  the  locking  limits.  When  A  is  small  the  oscillator  is  locked.  Periodic  pulling  occurs  just  beyond 
the  locking  band  (here  the  relative  phase  between  the  modes  oscillates  in  time).  Finally,  for  larger  A, 
mode  beating  occurs  with  a  continual  increase  in  relative  phase. 

When  F  <  0,  the  locking  band  can  become  very  large  since  x®  becomes  large  due  to  the 
subservient  nature  of  mode  2.  As  the  detuning  is  increased  (at  a  fixed  value  of  F),  the  amplitude  of  the 
damped  mode  decreases,  and  remains  phase  locked  to  the  growing  mode.  This  is  understandable  since 
the  damped  mode  has  no  source  of  excitation  but  the  one  driving  it.  Hence  the  two  modes  have  related 
phases.  In  this  case  the  predictions  of  the  linear  theory  are  relevant. 

The  linear  theory  adequately  predicts  the  values  of  x  and  O  in  the  locked  state  for  parameters 
well  within  the  phase  locking  band.  A  comparison  between  the  predictions  of  X°  and  <I>°  with  and 
without  self-saturation  is  shown  in  Fig.  10.  This  figure  again  is  based  on  gyrotron  parameters  chosen  in 
the  region  where  the  TE-ju  mode  has  positive  growth.  The  linear  results  are  reproduced  from  Fig.  8  for 
convenience.  The  self-saturation  coefficients  are  chosen  to  both  equal  -  1.0  x  10*^.  As  indicated 
previously,  these  locking  results  are  not  sensitive  to  the  precise  value  of  p.  The  shading  in  the  figure 
shows  the  region  of  unlocked  oscillator  operation.  In  the  gyrotron  case  considered,  the  values  of  the 
linear  cross-coupling  coefficients  are  so  small  that  the  modes  unlock  almost  as  soon  as  F  >  0  (this 
happens  at  a  magnetic  field  of  1.847  kG).  It  is  clear  that  the  predictions  of  linear  theory  are  inadequate  in 
this  region  of  parameter  space.  However,  for  F  <  0,  the  linear  theory  is  adequate.  Both  the  predictions  of 
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amplitude  ratio  and  locked  phase  become  more  accurate  at  high  magnetic  field  strengtiis  (larger  values  of 

iri). 

An  example  of  the  temporal  evolution  of  the  amplitude  ratio  and  relative  phase  is  given  in  Fig. 
11.  The  parameters  in  this  case  correspond  to  gyrotron  operation  at  1.853  kG  =  1.49,  F  =  -1.32, 

=  -0.0213,  C2  =  -0.0.532,  A  =  -0.188);  a  magnetic  field  for  which  Fig.  10  indicates  that  the  linear 
prediction  of  the  steady  state  is  not  a  good  one.  The  P's  are  chosen  to  be  -1.0  x  10'^.  An  interesting 
feature  of  this  figure  is  that  there  are  two  plateaus  in  both  the  x  and  <I>  curves.  One  corresponds  to  the 
true  nonlinear  phase  locked  state  (for  T  >  25),  while  the  other  (5  <  T  <  15)  is  a  "metastable"  equilibrium 
which  is  due  to  the  linear  locking  (and  occurs  at  the  values  of  and  <1>®  given  by  the  linear  theory  in  Fig 
10).  In  this  instance  the  self-saturation  is  small  enough  that  the  modes  achieve  the  linear  locked  state  but 
are  disrupted  as  the  modal  powers  increase  to  such  a  point  that  the  self-saturation  beconnes  important. 
The  linear  locking  result  has  validity  if  the  oscillator  pulse  terminates  before  the  metastable  equilibrium  is 
destabilized.  In  this  case,  the  upper  limit  on  pulse  length  is  about  20  e-folding  times  for  the  TE|j|  mode 

or  13.5  ns. 

When  F  >  0,  the  locking  bandwidth  generally  increases  with  an  increase  in  I F— 1|.  This  result  is 
shown  in  Fig.  12  for  the  parameter  set  (C|  =  0.8,  C2  =  0.5,  Pi  =  P2  =  -0.001).  The  gyrotron  example 

does  not  clearly  exhibit  this  effect  since  the  Cs  are  so  small  that  the  locking  band  becomes  extremely 
small  for  F  >0  and  is  difficult  to  compute.  The  shading  in  the  figure  indicates  the  region  over  which  the 
two  modes  are  phase  locked.  The  locking  band  is  asymmetric  about  A  =  0  for  the  same  reasons  given  in 
the  linear  case.  Fig.  9(a).  The  numerically  computed  locking  limits  can  be  checked  by  a  comparison  with 
(12),  using  the  associated  values  of  x^.  Good  agreement  is  found. 

Minima  in  the  locking  band  edges  are  seen  in  Fig.  12  at  F  =  0  (A  >  0)  and  F  =  2  (A  <  0).  These 
minima  occur  when  x®  reaches  such  a  value  that  the  denominator  in  (12)  is  minimized.  In  such  a  case  the 
detutung  can  be  varied  over  only  a  small  range  before  the  inequality  in  (12)  is  violated.  It  can  easily  be 
shown  that  x®  should  equal  one  for  a  minimum  to  occur.  At  this  amplitude  ratio,  (12)  indicates  that 
I  A/C|  I  =1,  at  the  locking  limit.  It  is  interesting  that  only  the  coefficient  C]  plays  a  role  here.  Equation 
(13c)  reveals  that  at  a  minimum  of  the  locking  edge:  COsO®  =  -  A/C|.  Using  this,  it  is  found  that  a 
phase  locked  equilibrium  of  (10),  (with  P's  included),  satisfies: 


Pmin  = 


F-1±2C2 

Pi  -p2 


Pi  *  P2 


(19) 


at  the  locking  edge  minimum.  Here  P^iin  normalized  power  in  each  mode.  The  two  different 
signs  correspond  to  the  locking  edge  with  A  >  0  (plus  sign)  or  A  <  0  (minus  sign).  When  the  self- 
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saturation  coefficients  are  equal,  then  a  condition  is  obtained  for  F  at  which  the  locking  minimum  will 
occur 


rMiN  =  -(-l  ±  2  C2),  Pi  =  P2  •  (20) 

The  two  signs  have  the  same  meaning  as  in  (19).  Note  that  both  (19)  and  (20)  depend  only  on  the 
coefficient  C2.  Equation  (20)  accurately  predicts  the  minima  in  Fig.  12  with  C2  =  0.5.  The  points  at 
which  it  is  most  difficult  to  phase  lock  (these  minima)  occur  very  near  F  =  1,  when  the  cross  coupling  is 
small.  This  is  just  the  point  at  which  linear  theory  predicts  an  infinite  lock-in  time. 

From  the  previous  analysis,  it  is  expected  that  the  locking  width  increases  with  an  increase  in 
I  I .  A  change  in  does  not  affect  the  size  of  the  locking  minima,  but  does  change  the  values  of  F  for 
which  the  minima  occur.  Hence,  for  a  fixed  ratio  of  growth  rates  Fg,  the  locking  band  could  either 
decrease  or  increase  with  a  change  in  depending  on  whether  approaches  or  departs  from  Fg. 

An  increase  in  the  absolute  size  of  the  self-saturation  has  no  effect  on  the  phase  locking  band  in  the  case 
Pi  ~  P2'  **  because  x°  depends  on  the  ratio  of  P|  to  P2  .  Unequal  saturation  coefficients  may  cause 

a  substantial  change  in  the  locking  bandwidth.  For  example,  if  <  1  ,  then  an  increase  of  Pi  causes 
mode  1  to  saturate  more  quickly  (and  at  a  lower  power  level)  so  that  x°  decreases  even  further.  This 
increases  the  denominator  of  (12)  allowing  a  wider  locking  bandwidth  (on  the  edge  where  x°  <  1).  Since 
it  may  not  be  true  that  x°  <  1  on  both  locking  edges  (A  <  0  and  A  >  0)  the  overall  locking  bandwidth 
measured  at  constant  F  may  either  increa^  or  decrease.  An  increase  of  P2  causes  the  opposite  effect. 
The  same  changes  in  Pi,  and  P2  have  the  opposite  effect  on  the  edge  of  the  locking  band  where  x°  >  1. 

It  can  be  shown  that  when  the  steady  state  amplitude  ratio  becomes  either  very  large  or  very 
small  the  locking  limits  become  infinite.  When  x®  »1,  (10)  yields  approximate  results  for  the  modal 
power  levels: 


Pi  =  -  —  and  P2  = - - — 

Pi  4A^Pi 

and  the  locking  criterion  (12)  becomes 


VCi2  +  €2^ 


(21) 
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which  is  clearly  always  satisfied.  For  x®  « 1,  the  power  levels  are 


P2  =  -  ^  and 

P2 


rci^ 

4A^  p2 


which  result  in  the  same  locking  criterion  as  (21).  An  example  of  this  phenomenon  is  shown  in  Fig.  12 
near  F  =  0.  It  can  be  seen  that  the  locking  edge  for  A  <  0  approaches  negative  infinity.  It  is  obvious  from 
the  figure  that  the  locking  edges  for  positive  and  negative  A  need  not  — >  for  the  same  value  of  F. 
Figure  12  shows  an  increase  in  the  locking  band  at  large  F.  This  is  due  to  a  decrease  in  x®  with 
increasing  F  at  both  edges  of  the  locking  band.  The  same  result  is  obtained  for  F  «1. 


V.  Conclusion 


The  system  of  two  electromagnetic  modes  interacting  through  amplitudes  and  phases  in  the  first 
order  approximation  has  been  analyzed.  Regimes  of  phase  locking,  periodic  oscillation  of  power  and 
phase,  and  mode  excitation  have  been  identified.  The  approximate  time  for  the  phase  locking  process  has 
been  found.  In  the  case  of  small  nonzero  (negative)  self-saturation,  it  is  found  that  the  phase  locking 
bandwidth  is  limited  in  size.  The  ratio  of  the  equilibrium  power  levels,  the  lock-in  time,  and  the  lock-in 
phase  are  fairly  well  approximated  by  the  linear  theory  near  the  center  of  the  locking  band.  C^alitative 
features  of  the  phase  locking  band  can  be  analytically  computed  and  have  been  numerically  verified. 
Further  work  is  proceeding  to  further  study  the  case  of  unequal  degrees  of  modal  self-saturation 

(Pi  ^  P2)- 
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Figure  Captions 


1:  Real  part  (dash)  and  imaginary  part  (solid)  of  (a)  Linear  cross-coupling  and  (b)  linear  growth 
coefficients  for  TE|n  mode  in  competition  with  TE||2-  **  the  detuning  of  the  TEm  frequency 

from  the  cyclotron  frequency. 

2:  Real  part  (dash)  and  imaginary  part  (solid)  of  (a)  Linear  cross-coupling  and  (b)  linear  growth 
coefficients  for  TE||2  mode  in  competition  with  TE^.  A(.2  is  the  detuning  of  the  TEn2  frequency 

from  the  cyclotron  frequency. 

3;  (a)  Imaginary  and  (b)  real  parts  of  cross-coupling  coefficients  of  TEjn  (solid)  and  TE||2  (dash) 
modes  as  a  function  of  the  static  magnetic  field. 

4:  Cross-coupling  between  the  TEjn  mode  and  (a)  TE113  and  (b)  TE114  modes.  Real  (dash)  and 

imaginary  (solid)  parts  are  shown. 

5:  Degree  of  cross<oupling  of  TEm  with  higher  order  modes. 

6:  Time  required  for  phase  locking  to  be  achieved.  Both  numerical  (solid)  and  exponential  decay  (dash) 
solutions  are  shown. 

7:  Imaginary  part  of  linear  growth  rate  of  TEm  (solid)  and  TE112  (dash)  modes  as  a  function  of 
magnetic  field. 

8:  Stable  phase  locked  amplitude  ratio  and  locking  angle  for  gyrotron  with  positive  TE|||  growth  rate. 

9:  Dependence  of  time  taken  to  phase  lock  on  (a)  mode  frequency  detuning,  (b)  relative  growth  rate,  (c) 
linear  cross-coupling. 

10:  Comparison  between  linear  theory  (solid  and  dash)  and  nonlinear  theory  (chaindot  and  dot)  with 
=  P2  =  -.0001. 

11:  Case  of  small  self-saturation  where  a  metastable  equilibrium  is  found  (p|  =  P2  =  --OfXlOl). 

12:  Phase  locking  bandwidth  in  the  nonlinear  case  as  a  function  of  relative  growth  rate  (P|  =  P2  “  -.001). 
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PART  II 


Coupling  Between  Low  Q  Modes  in  a  Beam-filled  Resonator 

with  Apertures 

Abstract  -  Part  II 

Mode  conversion  and  coupling  are  serious  problems  in  low  Q  gyrotron  operation,  as  they  can  result  in 
efficiency  reduction,  spectral  degradation,  and  undesirable  thermal  loading.  The  proper  modeling  of  this 
phenomenon  is  thus  an  important  goal  in  gyrotron  theory.  In  this  work,  the  electrodyrumics  of  a  low  Q 
gyrotron  cavity  is  described  in  terms  of  its  short  circuit  eigenmodes.  The  correct  set  of  eigenmodes  in  a 
simply  connected  gyrotron  cavity  are  identified  and  computed.  The  oscillator  equations  describing  the 
evolution  in  time  of  these  eigenmodes  are  derived.  Power  generation  and  loss  mechanisms  in  the  cavity- 
waveguide  structure  are  discussed.  The  electromagnetic  solution  is  verified  for  the  case  of  a  rectangular 
waveguide  coupling  out  of  a  cavity  of  rectangular  cross-section.  The  case  of  a  single  frequency  wave, 
propagating  along  the  waveguide  and  exciting  the  cavity  is  discussed,  and  is  compared  with  analytic 
theory  in  the  limit  of  a  short  circuited  waveguide.  The  theory  of  electromagnetic  waves  inside  a  cavity 
with  Ohmic  wall  losses  and  external  losses  is  investigated.  Implementation  of  these  developments  into 
multimode  gyrotron  theory  is  discussed. 
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1.  INTRODUCTION 


In  the  ensuing  work,  the  electromagnetic  fields  in  a  cavity  are  expanded  in  terms  of  the  short 
circuit  cavity  eigenmodes.  The  short  circuit  modes  correspond  to  the  cavity  completely  enclosed  by 
perfectly  conducting  walls.  These  short  circuit  cavity  eigenmodes  are  orthonormal,  and  form  a  complete 
set  with  respect  to  the  expansion  of  any  arbitrary  vector  field  in  the  cavity  volume.  Then  the  problem  of 
representing  the  fields  in  the  cavity  reduces  to  determining  the  expansion  coefficients  multiplying  the 
cavity  eigenmodes. 

The  problem  of  electromagnetic  field  expansion  in  a  cavity  not  completely  enclosed  by  a  perfect 
conductor  has  been  tackled  by  many  workers  over  the  years.  The  case  of  a  cavity  completely  enclosed  by 
perfectly  conducting  walls  was  solved  by  Weyl.^  Slater^  gave  a  comprehensive  theory  covering  this 
subject,  but  the  set  of  vectors  used  by  him  are  not  correct  in  a  cavity  with  open  boundaries;  he  imposed 
the  wrong  set  of  boundary  conditions  on  his  irrotational  functions,  and  as  a  result  missed  a  capacitive 
term  in  his  expression  for  input  imp>edance.  Moreover  his  set  of  vectors  correspond  to  mixed  boundary 
conditions  and  are  difficult  to  compute.  Teichmann  and  Wigner^  showed  for  lossless  cavities  that  an 
irrotational  magnetic  field  is  needed  to  completely  expiand  the  electromagnetic  fields  in  the  cavity  exdted 
through  holes.  This  irrotational  magnetic  field  vector  is  represented  as  the  gradient  of  a  scalar  function  u. 
The  function  u  is  determined  by  the  value  of  the  normal  component  of  the  magnetic  field  in  the  aperture, 
and  by  using  the  Neumann  function  of  the  cavity.  The  theory  of  Neumann  functions  and  their  properties 
is  covered  in  some  detail  by  Bergman  and  Schiffer,^  and  O.D.Kellogg.^  As  is  the  case  with  most  Green's 
functions,  the  Neumann  function  for  the  cavity  with  a  hole  in  the  surface  is  difficult  to  evaluate,  and  is  the 
main  drawback  of  the  theory  presented  by  Teichmann  and  Wigner. 

Schelkunoff^  has  demonstrated  that  Teichmann  and  Wigner’s  comments  about  the  short  circuit 
modes  being  incomplete  with  respject  to  the  expansion  of  electromagnetic  fields  were  not  correct  if  use  is 
made  of  a  short  circuit  which  conforms  to  the  impressed  field.  Borgnis  and  Papas^  wrote  a  very 
comprehensive  review  article  on  electromagnetic  waveguides  and  resonators,  in  which  they  present  the 
theory  of  cavity  field  expansions,  and  give  a  formula  for  computing  the  input  admittance  of  the  cavity 
resonator  at  the  cavity-waveguide  terminal  plane.  The  expansion  of  cavity  electromagnetic  fields  is  also 
treated  by  Kurokawa,®  and  indeed  this  p>ap>er  gives  the  most  correct  set  of  expansion  vectors. 
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II.  BACKGROUND 


In  this  section  some  of  the  main  ideas  in  the  theory  of  cavity  field  expansion  are  discussed.  First, 
the  complete  set  of  vector  fields  required  to  expand  the  electric  and  magnetic  fields  in  a  cavity  with 
apertures  on  its  surface  is  discussed.  Using  Maxwell's  equations  and  the  field  expansions,  coupled 
expressioirs  for  the  field  expansion  coefficients  are  derived.  Next,  the  oscillator  equations  which  describe 
the  time  evolution  of  each  of  the  cavity  modes  are  derived. 


A.  Complete  Sets  ofOrthonorma.  Vectors 

Kurokawa  has  shown  that  any  arbitrary  scalar  function  can  be  expanded  in  terms  of  an 
appropriate  set  of  complete  and  orthonormal  eigenfunctions  in  a  volume  V,  enclosed  by  a  surface  S,  if  the 
set  of  eigenfunctions  is  composed  of  the  solutions  of  the  wave  equation  with  the  Dirichlet  boundary 
condition: 


V2v„+k2  =  0  in  volume  V,  Va  *  ^  o*'  surface  S. 

The  set  Va  is  a  complete  set  of  orthonormal  functions  (a=l,2,3...),  and  can  be  used  in  expanding 
an  arbitrary  piecewise  continuous  scalar  function  in  V.  Another  complete  set  of  orthonormal  functions  of 
equal  generality  of  application  is  given  by  the  solutions  of 

V^<|iB+k3  =  0inV,  ^  =  0onS 
^  ^  9n 

where  the  derivative  is  taken  in  a  direction  normal  to  the  surface  S. 

An  arbitrary  piecewise  continuous  vector  function  in  the  volume  V,  can  be  expanded  in  a  similar 
manner.  One  set  of  complete  orthonormal  vector  functions  is  obtained  as  solutions  of 

V2ipp+k2?p  =  0  in  V.  nx  ?p=0.  V-  ?p=0  on  S  (p=0,l,2,...).  (1) 

-♦ 

Here  n  is  a  unit  outward  normal  from  the  cavity  surface.  Another  complete  set  of  eigenvectors  is 
obtained  as  solutions  of 

V2<l>q+k2«i)q  =  0  in  V.  iTx  Vx  <l>q=0,  n-  Oq=0  on  S  (q=0,lZ.  ).  (2) 
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Some  of  the  'Fp  and  have  the  same  eigenvalues  these  are  related  to  each  other  through  the 
curl  operation  as 

Vx'P()[=k(j<t>G[  and  • 

All  the  other  'Fpand  <l>q  that  do  not  have  common  eigenvalues  are  irrotational;  if  they  are 
denoted  as  and  G;^^(a=0,l,2..)  and  (X=0,1^...)/  then 

VxF£j=0  and  Vx  G;^=0. 

If  Fq  and  G;;^  have  nonzero  eigenvalues,  then  they  can  be  represented  through  the  gradient 
operation  as 


kaFa=Vva  and  k^  Gx=V<t)x 

where 

V^(j+ka  =  0  in  V,  Va  =  0  on  S 

and 

V^<|)X+kx  <t>X=0  in  V  and  ^-^=0  on  S  . 

dn 

The  vector  electric  and  magnetic  fields  in  the  cavity  are  expanded  in  a  manner  analogous  to  the 
above  expansions.  The  solenoidal  vectors  required  in  the  electric  field  expansion  are  denoted  as  E(x  and 
the  irrotational  vectors  in  the  electric  field  expansion  are  denoted  as  F^^  .  In  a  simply  connected  cavity 

however,  no  F^j  are  required  to  expand  the  electric  field  since  all  the  F„  are  orthogonal  to  the  real 

— ♦ 

electric  field.  That  this  is  so  can  be  quickly  seen  by  integrating  the  scalar  product  of  total  electric  field  E 
with  a  particular  modal  vector  Fjj  over  the  cavity  volume.  Assuming  that  the  actual  electromagnetic  field 
is  periodic  with  a  finite  period  CO,  this  leads  to  the  surface  integral  of  nxF„-H,  and  this  vanishes  because 
from  eqn.(l),  it  can  be  seen  that  nxFjj  is  equal  to  zero  on  the  surface  of  the  cavity. 

In  a  similar  fashion,  the  magnetic  field  is  expanded  in  terms  of  the  solenoidal  vectors  and  the 
irrotational  vectors  Gx  •  That  the  irrotational  vectors  Gx  are  required  to  expand  the  magnetic  field  may 
look  strange,  but  that  this  is  the  case  can  be  seen  from  the  Helmholtz  theorem,  which  states  that  a  given 
vector  field  in  a  region  V  enclosed  by  a  surface  S  can  be  expressed  as  the  sum  of  the  gradient  of  a  scalar 
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field  and  the  curl  of  a  vector  field.  Only  when  the  divergence  of  the  given  vector  field  is  zero  in  V,  and  its 
normal  component  zero  on  S,  can  it  be  set  equal  to  the  rotation  of  another  vector  field.  Similarly,  only 
when  the  curl  of  the  given  vector  field  is  zero  in  volume  V,  and  its  tangential  component  equals  zero  on 
the  surface,  can  it  be  set  equal  to  the  gradient  of  a  scalar  field.  In  the  first  case,  the  vector  field  is  purely 
solenoidal,  and  in  the  second  case  it  is  purely  irrotational.  The  boundary  condition^  on  the  surface  are 
important  to  distinguish  the  purely  solenoidal  or  purely  irrotational  vector  fields,  and  if  a  vector  field  is 
not  purely  irrotational,  then  even  if  its  curl  is  zero  in  the  volume  V,  it  can  not  be  generally  written  as  the 
gradient  of  a  scalar  field.  In  the  same  way,  even  if  the  divergence  of  a  vector  field  is  zero  in  the  volume  V, 
if  it  is  not  purely  solenoidal,  then  it  can  not  generally  be  written  as  the  curl  of  a  vector  field.  Thus  we  can 
write 


E  =  IEa  EE„dv 
a  Jv 


(3) 


H  =  IH„ 
a 


HH^dv+XGx  HGxdv 
V  ^  Jv 


(4) 


3=IJ"‘ai  JE(xdv 
a  Jv 


(5) 


P=  IVa 
a  Jv 


(6) 


where  .1  and  p  are  the  volume  current  and  charge  densities  respectively.  All  the  integrations  in  eqs.(3)- 
(6)  are  done  over  the  volume  of  the  cavity.  The  eigenmodes  satisfy  the  relations 


V^Ea+kJ  E„  =  0  in  V,  n  X  Eot=0.  V-  £^=0  on  S+S' 

(7) 

V^H{3(+k„  Hq[  =  0  in  V,  iTx  Vx  Hq[=0,  n-  Hq(=0  on  S+S 

(8) 

V^Gx+k^Gx  =  0  in  V,  iTx  Vx  Gx*0.  n  Gx=0  on  S+S' 

(9) 

VxEQj=lq3[H(3[  and  VxHQj=k(xEQ[  and  VxG^  =  0 

(10) 

and  the  orthonormality  relations 
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(11) 


j  J  J  ^a’®P*^~^ap 

where  6(xp  is  the  Kronecker  delta  function  (6o(p=l  if  a  =  P  and  5(xp=0  if  a;*  P).  Here  S  denotes  the  portion 
of  the  cavity  surface  which  is  a  conducting  nraterial,  and  S'  denotes  the  portion  that  is  an  aperture. 

From  Maxwell's  equations  it  is  seen  that  VxE  has  the  properties  of  a  magnetic  field,  hence 
and  are  used  in  the  expansion  of  VxE.  Similarly,  since  VxH  has  the  properties  of  an  electric  field, 
E(2  are  used  in  its  expansion.  Thus 

VxE  =  Z  Hq(  VxE-  H^dv  -«■  £  j  VxE-  G^dv 

a  Jv  ^  Jv 

and 

VxH  =  I  Vxii- E„dv. 
a  Jv 


B.  Oscillator  Equations 

From  vector  identities  and  using  the  properties  of  the  eigenvectors,  the  following  relations  are 
obtained 


VxE-  Ha  =  V.^  Ha  kaE-  Ea 
VxE-  G;^  =  V-Ex  G;;^ 

VxH-  Ea  =  V-Hx  Ea  +  kaH-  Ha  . 

Substituting  these  relations  in  eqs.  (3)  &  (4),  using  the  divergence  theorem ,  eqn.(ll),  and  using 
the  boundary  conditions 

nxEa=0  on  S■^S'  andiT  Ha=0onS-i-S' 


the  following  equations  are  obtained 


k(x|  E-E(jdv  +  ji— I  H-  Ho[dv  =  -j  nxE-  H(jC 

Jv  Vv  Js.s’ 

fil  I  H.  G^dv  +  I  nxE  Gxds  =  0 . 

•V  »s+s 


The  following  two  oscillator  equations  are  obtained  by  eliminating  the  expansion  coefficients 
from  the  coupled  eqs.  (12)-{14) 


-L-^f  E  Ettdv  +  k^  EiLdv*-k„  iixEHttdS-p— [  J  Eadv  (15) 

c^dt^lv  Jv  is+s' 

-  HHojdv  +  kJ  HHadv  =  kci  JEadV-e—  iixHHads  (16) 

Jv  Jv  Jv  ^JS+S’ 


J-^f  H  H„dv  +  k2f  H  Hadv  =  k„f  JEadV-e— f  iixH  Hads  (16) 

c2at2jv  iv  “  “iv  dJs+s- 

These  oscillator  equations  form  the  basis  of  the  subsequent  gyrotron  cavity  theory.  They  will 
now  be  solved  subject  to  different  boundary  conditions. 


m.  EXQTATION  OF  CAVITY  FIELDS 

In  this  section,  the  oscillator  eqs.(15)  and  (16)  are  solved,  keeping  different  source  terms  on  the 
right  hand  side  in  each  problem.  This  is  a  useful  exercise  because  it  enables  us  to  check  the  results  of  the 
field  expansion  theory  in  cases  with  relatively  simple  geometry.  By  including  more  and  more  terms  on 
the  right  hand  side,  we  make  an  orderly  progress  toward  the  gyrotron  problem,  which  is  solved  using  a 
circular-cylindrical  geometry  and  a  convection  current  in  the  interior  of  the  cavity. 

A.  Cavity  Excited  by  a  Waveguide 

The  oscillator  eqs.(15)  and  (16)  together  with  eqs.(3)  and  (4)  describe  completely  the  time 
evolution  of  electromagnetic  fields  in  a  cavity.  Here  as  a  first  step,  these  equations  are  solved  in  a  cavity 
excited  by  a  waveguide  (Fig.l)  when  a  single  frequency  wave  propagates  down  the  waveguide.  This 
steady  state  problem  is  solved  in  the  frequency  domain  by  performing  the  Fourier  transformation  of  the 
oscillator  equations .  Writing 
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f(X=j  E-  EqcIv  and  h{x=|  H-  H^dv  and  g;^^=|  H-  G;^^dv 

Jv  Jv  Jv 


and  Fourier  transfoiming  eqs.(12),  (14  )  and  (15),  one  can  write 


M“)  =  '*^a(“)  and  gx(w) 
ka 


nxE(r.(o)-  G;^_ds 


-  ^jfa{“)  =  ->^1^  ^  nxE(r.0))-Ha(r)d5  +  j coji  J  jfr.o))-  Ea{r)dv 


When  it  is  assumed  that  there  is  no  beam  in  the  cavity,  then  the  second  term  on  the  right  hand 
side  vanishes.  If  it  is  also  assumed  that  the  walls  of  the  cavity  are  made  of  a  perfect  conductor,  then  the 
tangential  component  of  the  electric  field  is  equal  to  zero  on  S,  and  eqn.(18)  becomes 


|4  -  ^Jfa{“)  =  j  nxE{r.(o).  hJ 


To  solve  for  f„((o)  we  need  a  knowledge  of  iwE  at  the  waveguide-cavity  interface  (aperture). 
This  field  may  be  found  by  matching  the  field  in  the  waveguide  with  that  of  the  cavity  across  the 
aperture.  The  waveguide  fields  expressed  as  summation  of  waveguide  eigennxxles  are 


H  =  IisGt.s+^z,s 

S  4u»C 


where  s  represents  the  s'th  waveguide  mode,  and  t  and  z  denote  the  transverse  and  the  longitudinal  (in 
the  increasing  direction  along  z-axis)  directions.  F(j5  and  G  satisfy  the  orthonormality  relations  given  by 


j  Ftr-  F,sds=j  Gtr  Gtsds=|  FyxGts  z"ds  =  6re 

Js'  Js'  Js' 


where  6^,  is  the  Kronecker  delta  function  as  defined  earlier,  and  z  is  a  unit  vector  in  the  z-direction. 
Substituting  eqs.  (20)  and  (21)  into  eqn.(19)  results  in 
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where 


fa((0)  = 


-^lYosVs 

1.2  0)2  s 
“*T 


Jos 


-  I  ^a‘ 
Js' 


(23) 


For  the  rectangular  cavity  shown  inhg.l,  the  TE(to  z)  cavity  eigenmodes  are  as  follows: 


H(j,  =  y  H(v  Ert  =  -ja)(3il^^  , E(y  = 

^  az2  ^  axaz  “y  ayaz  “*  “y  “^T 


axaz  ’  “y  "  ayaz 


where 


a  denotes  the  triple  indices  (m,n,p)  and  €  and  £  p  are  the  Neumann  factors  which  take  on  a  value 
of  2  if  the  subscript-index  is  not  equal  to  zero,  and  have  a  value  of  1  if  the  subscript-index  equals  zero. 
Similarly,  the  TM(to  z)  cavity  eigenmodes  are  dehned  as: 


where 


Va  = 


abd 


E 


with  h,  k(3[ ,  and  o)(x  defined  as  above.  The  G^,  are  defined  as 


Gi  =-L^,Gi  =-L^,Gi 
^  kx  ax  T'  Irn  ;)v  ^ 


-  1 
kx  ay 


kx  dz 


where 
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and 


The  waveguide  eigenmodes  are  defined  in  a  similar  way.  The  TE(to  z)  eigenmodes  are  defined 


Eg  =  -jconVx  Yj,  where  Vs  =  e'fP^ 


Hg  =  VxVx  Vs.  where  Vs  =  cos|^^j  e'jP^ 

where  s  denotes  the  double  waveguide  indices  (m,l),  e  e  j  are  the  Neununn  factors  as  defined  earlier, 

h  and  ^  are  defined 

h  =  and  P  = 


The  TM(to  z)  eigenmodes  are  defined 


E  =  Vx Vxv,  where  v  =  sin^Da^J sinj^^j e'iP^ 


H  =  j  (xV^,  where  v  =  -z.  ^  ^  sin^ffl^j  sin|^^j  e'JP^. 

The  Vj  are  the  amplitude  coefficients  of  waveguide  electric  field  modes  at  the  aperture;  they  need 
to  be  computed  before  fQ((a)  can  be  evaluated.  From  eqs.(20)  and  (21)  it  can  be  seen  that 

ir=  HG^ds.  (24) 

Js' 


V  c2  a2 


By  substitution  of  eqn.(5)  into  eqn.(24) 

>r=S  YfgVs 
s 


where 


‘rs 


-ja>  V  'ygsTm- 
^  ®  (0^-a)2 


+ 


r^IPXfPXs 

J«fi  A, 


(25) 

(26) 


44 


with 


P3lr=f  Gx-Gyds. 

Js' 

The  Yn;  are  admittance  coefficients,  and  describe  how  the  current  in  the  r'th  mode  of  the 
waveguide  couples  into  various  cavity  modes,  and  how  each  of  them  in  turn  couple  back  into  the  s'th 
waveguide  mode.  The  above  equation  for  if  actually  describes  a  system  of  linear  equations  which  are 
solved  to  compute  the  Vj..  The  solution  is  given  in  the  following. 

A  small  distance  away  from  the  aperture,  there  are  a  great  number  of  waveguide  modes 
propagating  in  both  forward  and  backward  directions,  because  of  the  impedance  mismatch  at  the 
aperture.  While  the  amplitude  of  most  of  these  modes  will  fall  off  rapidly  away  from  the  aperture 
(because  the  waveguide  is  designed  for  dominant  mode  propagation),  their  amplitudes  right  at  the 
aperture  are  not  necessarily  negligible.  Thus  the  Vj.  at  a  point  in  the  waveguide  is  composed  of  v^  and 
with  the  superscripts  indicating  waves  propagating  in  the  +z  and  -z  directions. 

Vj5  =  v^-V5  andij;  =  i^  +  ij5*  for  TE  modes  (27) 

and 

Vjj  =  v^  +  Vg  and  ijj  *  i_^  -  i^  for  TM  modes  (28) 

with  i^  =  Y(jj^Vs  and  ig  =  where  is  the  admittance  of  the  s'th  mode  . 

Using  matrix  notation,  eqn.(25)  can  be  written  as 

('^Te)  =  (YoarEj'^iYiE-TE)  ({''te)  '  (’'iE)HYTE-TM)  (('^Tm)  ('^Tm))]- 

and 


('^Tm)  ■  ('^Tm)  =  (Yo>rM )  YtM-Te)  {(^'Te)  '  (''te)M Ytm-tm)  ((^tm)  +  (''m))]- 

upon  rearrangement,  these  equations  yield 

(^te)  =  (^1 )  (^te)  (^2)  (vtm)  (''Tm) = (®1  ){vte)  (®2 )  (vtm)  • 

Thus  V jg  and  vjj^  are  solved  in  terms  of  the  known  quantities  vjg  and  v^.  From  eqn.(23),  it  can  be 
seen  that  all  the  f^  can  be  determined  when  the  v,  are  known. 
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The  results  obtained  by  solving  the  above  system  of  equations  may  be  best  checked  in  the  case  of 
a  short-circuited  waveguide  (Fig.  2)  by  making  the  height  of  the  cavity  and  height  of  the  waveguide 
equal.  It  is  assumed  that  the  waveguide  is  designed  for  the  dominant  TE^q  mode  propagation,  and  that 
the  amplitude  ^  shows  the  normalized  admittance  at  the  apjerture  as  a  function  of  the 

frequency.  The  result  obtained  from  the  cavity  field  expansion  theory  is  compared  with  an  analytic 
expression  for  admittance  at  the  aperture  plane  in  the  short-circuited  waveguide.  As  can  be  seen,  the 
results  are  in  very  good  agreement.  This  can  also  be  seen  by  comparing  the  expressions  for  ap)erture 
adnuttance  from  the  cavity  eigenmode  expansion  theory  and  the  short  circuited  waveguide. 

For  a  short  circuited  waveguide,  the  modal  admittance  for  the  donrunant  TEpo  mode  at  a  distance 
d  from  the  shorted  end  of  the  waveguide  is  given  as 

Y((d)  =  -j-^  cot{Pd)  where  P  =  'I  /  — —  —  (29) 

ft®  V  a^ 


Expanding  this  function  Y(a})  as  a  Mittag-Leffler  series  of  partial  fractions  using  the  theory  of 
residues,  we  obtain 


P=(  °  k^pid  0)^  - 


+  .-S-codiSd.. 
jcopa  a 


(30) 


Equation  30  is  the  same  as  eqn.(26),  when  eqn.(26)  is  expanded  over  the  TEpOp  typ^e  modes  (since 
in  the  shorted  waveguide  only  the  TEpop  cavity  eigenmodes  are  needed  in  the  expansion  for  the 
dominant  TEpo  waveguide  mode  pro{>agation).  The  Erst  sununation  on  the  right  hand  side  in  eqn.(30)  is 
the  contribution  of  the  solenoidal  modes  to  the  admittance  function,  and  the  second  term  is  the 
contribution  of  the  irrotational  modes  to  the  admittance.  (The  ap>erture  admittance  exhibits  resonance 
behavior  at  the  cavity  (lOp  modes)  resoiuint  frequencies.  This  can  be  readily  seen  from  eqn.(29)  and 
eqn.(30). 

Figures  4  and  5  show  the  y-compx)nent  of  the  electric  field  plotted  vs  z-axis  along  the  center  of 
the  waveguide  for  the  short-circuited  waveguide  case.  Again,  the  results  from  the  cavity  held  exp>ansion 
theory  are  compared  with  the  waveguide  theory  expression  for  electric  field  in  the  shorted  waveguide.  As 
can  be  seen,  the  results  are  in  very  good  agreement  except  very  near  the  apjerture,  where  the  field 
obtained  from  the  cavity  field  exp^ansion  drops  off  to  zero,  unlike  the  waveguide  field.  This  is  to  be 
expjected,  since  all  the  eigenmodes  used  in  the  cavity  field  exp»nsion  have  a  zero  tangential  component 
on  the  cavity  surface.  This  is  a  limitation  im|X)sed  by  the  nature  of  the  modes  used  in  the  exp>ansion,  and 
can  not  be  corrected  regardless  of  how  many  modes  are  used  in  the  exp^ansion.  However,  by  irKreasing 
the  number  of  eigenmodes  used  in  the  expansion  of  cavity  fields,  the  exp^anded  field  can  be  made  to 
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better  agree  with  the  true  field.  This  result  is  somewhat  analogous  to  the  Fourier  series  expansion  using 
sine  terms  only. 

Figure  6  shows  the  normalized  admittance  of  a  cavity  as  a  function  of  frequency  when  the  cavity 
and  the  waveguide  have  different  heights.  The  admittance  function  exhibits  resonance  behavior  at  all  the 
0)q(  (a  =  Inp  modes).  Since  the  width  of  the  waveguide  and  cavity  are  equal,  only  (a  =  Inp )  modes  need 
to  be  considered.  Figure  7  shows  the  profile  of  the  y-component  of  the  electric  field  along  the  cavity  axis, 
when  the  cavity  and  waveguide  have  different  heights.  As  can  be  expected,  the  field  drops  to  zero  near 
the  aperture. 

Figure  8  shows  the  energy  stored  in  the  cavity  as  a  function  of  the  frequency.  The  stored  electric 
energy  in  a  cavity  can  be  computed  from 


The  stored  magnetic  field  energy  is  computed  from 


f--  u 

• 

1  HHdvs(i 

+  gaG(x|- 

jlhpHp  +  gpGp 

Iv  ^ . 

Wm  =  ^  l(|ha|2+|gd2) 

4  a 

and  the  total  energy  stored  in  the  cavity  is  given  as 


W  =  We  +  W„,. 

For  the  purpose  of  comparison,  the  electric  stored  energy,  magnetic  stored  energy,  and  the  total 
energy  are  computed  for  the  short  circuited  waveguide  fields.  Figure  8  shows  the  total  energy  (computed 
from  both  the  cavity  expansion  fields  and  the  shorted  waveguide  fields)  as  a  function  of  frequency.  As 
can  be  seen  from  the  figure,  they  are  in  very  good  agreement.  The  resonant  frequencies  for  the  shorted 
waveguide  (equivalent  of  a  cavity  with  one  side  completely  removed)  are  shifted  down  from  the 
completely  enclosed  cavity  case,  and  this  is  as  expected  from  theory.  The  short  circuited  waveguide  is  a 
very  low  Q  cavity,  and  from  Fig.8  it  can  be  seen  that  all  the  Q  factors  have  magnitudes  of  the  order  of  10. 

B.  Cavity  with  Walls  of  Finite  Conductivity 

In  a  cavity  enclosed  by  walls  of  finite  conductivity,  the  field  profile  in  the  cavity  is  a  little  different 
from  the  case  of  a  cavity  made  of  a  perfect  conductor,  because  different  cavity  modes  are  coupled  through 
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the  surface  impedance  of  the  cavity  walls.  In  the  analysis  presented  here,  it  is  assumed  that  the  cavity 
wall  conductivity  while  finite,  is  still  very  large.  Using  this  assumption,  a  perturbation  analysis  is  done  to 
solve  the  problem. 

Writing  nxE(r.(a )  =  ti((o  )H(r,(o )  on  the  conducting  portion  of  the  cavity,  and  defining  the  Q  factors 

where  6  is  the  skin-depth  of  the  conducting  walls  given  by 

S=  y,/  — 2 —  and  (0®P  =  VWftWg  . 

V  (O^Ppa 


The  subscript  "o"  indicates  an  Ohmic-Quality  factor.  Thus  we  have 


Using  matrix  notation,  we  can  write 


(f)  =  (A)(f)+(B)(v). 

In  the  above  expression,  all  the  terms  that  have  a  ^  dependence  are  deemed  to  be  1st  order  terms 

(since  -1-  «  1)  and  terms  with  -L  dependence  are  the  second  order  terms  and  so  on;  otherwise  they  are 

Q  q2 

of  zeroth  order  (~  1).  Similarly, 
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+ 


Using  nutrix  notation,  we  can  write 


(i)  =  (C)(f)+(D)(v) 


An  ordering  of  the  teims  is  used  so  that 


f  =  +  5f,  V  =  yO  +  Sv,  i  =  i®  +  8i, 

A  =  A^  +  A^,  B  =  +  C  =  cP  +  cKand  D  =  dP  +  D^ 

and 

Sf-«  1,  5sl_«  1,  5L«  i,etc. 

fO  yO  iO 

Then  we  can  write 

(fi)  +  (6f>|(A0)f(A^ll  [(t®M8f)M(BO)f(B^)I(vO)+{8v)]  (31) 

and 

{i0)+  (Si>((cOMc^)]  [(^k8f)M(D°MD^)I(vOM8v)].  (32) 

Thus  we  have  2  sets  of  linear  equations  for  the  2  sets  of  unknown  variables  (5f)  and  (5v).  All  the 
perturbed  coefficient  matrices  (a^),  (b^),  (c*),  (d^)  have  teims  with  a  ^  dependence.  (5i)  is  linearly 

related  to  (Sv)  via  eqs.(27)  and  (28).  Equating  the  zeroth  and  the  first  order  quantities  separately  in  the 
above  system  of  equations,  the  quantities  (5v)  and  (5f)  are  found  in  terms  of  (f^)  and  (vB).  The  solution  is 
obtained  as 


{8f)  =  {Ei){fO)+(E2)(vO)  and  (8v)  =  {Fi)(f<^)+ (F2){v0). 

Here  (Ei),  (E2),  (Fi),  and  (F2)  are  coefficient  matrices  which  are  obtained  by  rearranging  eqs.(31)  and  (32). 
The  full  expressions  for  these  terms  are  sufficiently  complicated  as  to  preclude  their  inclusion  in  this 
report. 
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IV.  CONCLUSIONS 


The  correct  set  of  eigenvectors,  required  in  the  expansion  of  any  arbitrary  electronnagnetic  fields 
in  a  cavity  with  apertures  in  its  surface,  are  identified  and  computed.  The  oscillator  equations  that 
describe  the  time  evolution  of  the  amplitude  coefficients  of  the  eigenmodes  are  derived.  Using  the 
geometry  of  a  rectangular  cavity  fed  by  a  waveguide,  these  coefficients  are  computed  in  the  frequency 
domain.  A  short  circuited  waveguide  geometry  is  used  to  check  the  results  obtained.  The  theory  of  a 
cavity  enclosed  by  walls  of  finite  conductivity  is  investigated.  By  including  the  J-  E  term  on  the  right 
hand  side  of  the  oscillator  equations,  the  gyrotron  cavity  problem  can  also  be  solved. 
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HGURE  CAPTIONS 


1:  Geometry  of  the  waveguide  feeding  the  cavity  (dimensions:  a=0.01783m,  b=0.025m, 
c=0.00993m,d=0.04m). 

2:  Geometry  of  the  short  circuited  waveguide  (dimensions:  a=0.01783m,  b=c=0.(X)993m, 
d=0.04m). 

3.  Normalized  aperture  admittance  as  a  function  of  frequency  in  the  shorted  waveguide  for  the 
TEio  mode. 

4.  Electric  field  profile  along  the  cavity  axis  in  the  short  circuited  waveguide  for  the  TEiomode 
(waveguide-cavity  aperture  at  z=0,  shorted  end  at  z=0.04m). 

5.  Expanded  view  of  the  electric  field  profile  along  the  cavity  axis  in  the  short  circuited 
waveguide  for  the  TE^q  mode  (waveguide-cavity  aperture  at  z=0,  shorted  end  at  z=0.04m). 

6.  Normalized  aperture  admittance  of  the  TEjo  mode  as  a  function  of  frequency  for  the  cavity 
depicted  in  Fig.l. 

7.  Electric  field  profile  along  the  cavity  axis  for  the  TEio  mode  in  the  cavity  depicted  in  Fig.l . 

8.  Normalized  total  stored  energy  in  the  cavity  as  a  function  of  frequency  for  the  TEio  mode  in 
the  shorted  waveguide. 
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PART  III 


Experimental  Work  on  a  12  GHz,  TE^oi  -  TEgn 
Gyrotron  Oscillator 

Abstract  -  Part  III 

Results  and  discussion  are  given  on  the  design  and  test  of  a  two-mode  gyrotron  oscillator.  The 
predictions  of  the  quasi-linear  mode  coupling  theory  for  this  specific  case  are  given.  It  is  found  that  the 
coupling  will  be  quite  strong,  with  the  cross-coupling  coefficients  on  the  same  order  as  the  growth  rates. 
Substantial  modal  tuning  is  required  to  prevent  the  growth  rates  of  both  modes  from  being  identical.  In 
this  case  only  the  case  r=l  will  be  tested,  (this  allows  only  a  restricted  regime  of  predicted  phenomena  to 
be  examined).  The  resonator  design  is  outlined  and  cold  test  results  are  given.  HFSS  simulations  of  the 
three  dimensional  cavity  structure  and  coupling  probes  are  included.  It  has  been  found  that  the  coaxial 
probes  can  couple  very  strongly  to  the  cavity  (estimated  Q  on  the  order  of  100-200)  without  intercepting 
any  electron  beam.  The  design  and  calibration  of  beam  and  rf  diagnostics  is  also  described.  Electron 
beam  collector  design  and  the  design  of  iron  pole  pieces  for  the  solenoidal  magnet  are  given.  A  MAGIC 
simulation  of  the  propagation  of  high  and  low  energy  gyrating  electron  beams  through  the  collector 
region  is  described.  The  electron  gun  trim  coil  design  and  POISSON  simulations  of  the  trim  coil  and  the 
main  solenoid  are  given  and  compared  with  measurements. 
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I.  INTRODUCTION 


The  gyrotron  quasi-Iinear  theory  has  been  employed  to  describe  the  mode  competition  between 
two  closely  spaced  modes  in  a  rectangular  cavity.  A  code  has  been  written  and  growth  rates  and  cross- 
coupling  coefficients  can  be  computed.  A  sununary  of  results  of  this  investigation  are  given  in  Section  II. 
The  rectangular  cavity  and  tuning  structure  have  been  fabricated  and  cold  tested.  The  design 
considerations  for  this  cavity  are  included  in  Section  III.  Results  of  design  and  testing  of  the  experimental 
infrastructure  is  given  in  Section  IV.  This  includes  the  construction  of  a  beam  tube  (to  test  beam 
propagation  through  the  system),  electrostatic  beam  probes,  trim  coil  electromagnet,  and  electron  beam 
collector.  Photographs  of  these  devices  are  included. 


II.  RESULTS  OF  QUASI-UNEAR  THEORY 


The  computation  of  the  linear  cross-coupling  between  rectangular  TEiOq  and  TEoiq  modes  has 
been  carried  out  using  the  quasi-Iinear  approximatioiv  The  rate  equations  have  the  form: 


CD;  Pi 


Vctot 


Vj  =  “j  -  “i  -  2^  2  “i '+  2  (“i*  ®  -  “j* '*) 


Here  Pj  and  Vj  are  the  power  and  phase  of  the  )th  mode,  respectively.  Primes  denote  real  parts  end 
double  primes  denote  imaginary  parts.  The  phase  angle  -d  is  defined 

d=(cOj-(flk)t  +  Vj-Vk  • 

The  linear  growth  rate  and  frequency  pushing  aj  are  defined 


Yo 


1- 


Pl 


(YomcfJ  Tt^me, 


where: 
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The  coefficients  of  cross-coupling  are  defined 


tti  v  = 


1 


Yo 


Pi 


(Yo  m  c) 


8Ioe 


n^me. 


Jt‘0 


where: 


aj(4)=f  f  -f  e‘^^"fqd^”d^’ 

Ji'mO  Jf«0 


Figure  1  shows  the  linear  coefficients  for  the  TEioi  mode  of  the  cavity  to  be  described  later.  The  curves 
represent  a'  (solid),  a"  (dash),  a'loi^oil  (chaindot),  and  a"ioi  oil  (x’s).  Because  the  resonant  frequencies 
are  so  close  together,  the  coefficients  for  the  TEoii  mode  are  nearly  identical.  By  adding  a  mismatch  at 
the  cavity  output,  it  is  possible  to  increase  the  cavity  quality  factor,  (these  quality  factors  are  dominated 
by  external  loss).  This  variation  of  quality  factor  causes  the  coefficients  to  be  quite  different  for  the  two 
inodes.  It  is  of  interest  to  determine  the  range  of  validity  of  this  theory.  Hence  the  experiment  will  be 
compared  with  these  results.  More  accurate  nonlinear  computations  can  later  be  made,  and  a  further 
comparison  with  experiment  can  be  done. 


III.  GYROTRON  CAVITY  DESIGN  AND  COLD  TEST 


A.  Design  ConsideraHons 

The  design  of  the  two-mode  experiment  proceeded  as  follows.  Once  the  cavity  modes  are  fixed, 
there  are  two  free  parameters  in  our  design.  These  are  the  beam  energy  and  the  radiation  frequency.  All 
other  design  values  can  be  derived  from  these  two.  The  other  quantities  of  interest  are  the  magnetic 
field,  cavity  length,  beam  axial  velocity,  beam  a,  cutoff  frequency  of  waveguide,  beam  guiding  center 
radius,  cutoff  frequency  of  beam  tunnel. 

First  the  beam  energy  is  fixed.  The  hard  tube  modulator  in  our  lab  will  generate  beams  with  up 
to  40  keV  of  energy.  We  investigated  two  different  beam  energies  at  radiation  frequencies  between  12 
and  13  GHz.  The  frequency  range  is  chosen  because  of  the  ease  of  obtaining  components  and  the  ease  of 
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construction  of  the  cavity.  Four  equations  are  used  to  find  the  parameters  suitable  for  resonance  between 
the  beam  wave  and  electromagnetic  field.  These  are: 


0)  =  ^co2  +  k  (1) 

(ji)  =  Q/y  +  k^w^  (2) 


a  = 


1.5^677 


V^l  +  1.5^(l-B/7) 


(3) 


\^  =  csj  \-\hp-- (^f 


(4) 


The  first  two  equations  give  the  dispersion  relations  for  the  electromagnetic  wave  and  the  beam  wave 
while  the  third  equation  (written  for  the  magnetic  field  B  in  kG)  gives  the  beam  a  based  on  a  cathode 
magnetic  field  of  350  G,  cathode  radius  of  1565  cm  which  produces  a  beam  with  a  of  1.5.  Equation  (4) 
reflects  the  definition  of  y  in  terms  of  the  electron  velocities.  Qtmbining  eqns.  (3)  and  (4)  we  obtain  a 
relation  between  B  and  Further  requiring  that  the  wave  phase  and  group  velocities  be  matched  at  the 
resonance  point  results  in: 


v 


z 


2  +  k,2c2 


(5) 


Eliminating  oic  from  eqns.  (1)  and  (5)  we  obtain  v^  =  c2  Using  this  relation  to  elimiiuite  v^  from  our 
relation  between  B  and  Vz,  and  eliminating  kz  from  the  same  relation  by  using  vz  =  c2  kz/w  in  (1)  we 
arrive  at  a  solution  for  B: 

_ yf _ 

0.0989(1-1/72) 

where  r  =  (2  re)  2.8  x  10^  and  B  is  in  kG.  Now  the  other  parameters  can  be  found: 

\/ ^  =  ,  (0  =  v^Co2-k  ^c2 
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and  eqn.  (4)  for  vx  and  a  =  vx/vz.  From  this  analysis  it  is  found  that  for  high  voltages  (~  30  kV)  the 
cavity  length  is  quite  short  (~  5  cm),  the  a  is  low  (~  0.8),  and  the  beam  tunnels  are  not  well  cutoff  (fj.  ~  16.5 
GHz).  The  requirement  of  n\atched  phase  velocities  was  relaxed  (in  order  to  make  the  cavity  longer)  and 
it  was  found  that  competition  from  other  axial  modes  became  a  problem.  Short  cavities  have  the 
disadvantage  that  the  start  oscillation  current  can  be  quite  high  for  a  cavity  Q  of  ~500.  Thus  a  lower 
voltage  was  used.  For  a  beam  voltage  of  18  kV  and  radiation  frequency  12  GHz,  it  was  found  that  a  =  1.4, 
cavity  length  is  8.38  cm,  and  the  beam  tunnel  cutoff  is  19.87  GHz.  This  case  requires  a  magnetic  field  at 
the  cathode  of  225  G.  Figure  2  shows  the  phase  and  group  velocity  matching  and  the  separation  from  the 
next  higher  mode;  the  TE102. 

The  machine  drawing  for  the  final  cavity  design  is  shown  in  Figure  3  (all  dimensions  are  in 
inches).  The  cavity  proper  is  shown  at  the  bottom  of  the  figure,  with  nominal  dimensions  of  1.264  x  1.264 
cm.  The  cavity  length  is  8.38  cm  as  described  previously.  The  upper  side  of  the  cavity  can  be  moved  by  a 
tuning  plunger  to  allow  adjustment  of  the  cavity  resonant  frequency.  The  plunger  is  geared  to  the 
outside  by  an  anti-backlash  worm  gear,  a  worm,  and  a  025”  303  stainless  steel  shaft  which  feeds  through 
the  magnet  pole  piece.  The  gearing  is  such  that  movement  of  the  tuning  plunger  can  be  controlled  to 
within  0.10  pm.  This  accuracy  is  needed  since  the  cavity  resonant  frequency  is  a  sensitive  function  of 
cavity  transverse  dimensions.  The  cavity  cap  is  attached  via  twelve  number  6  screws  which  compress  an 
OFHC  rectangular  copper  gasket  of  80  mil  crossectional  diameter.  The  tuning  mechanism  is  supported 
by  an  aluminum  superstructure  attached  to  the  outside  of  the  cavity.  The  rf  seal  is  performed  by  a  2  mil 
thick  stainless  steel  wiper  which  is  electron  beam  welded  to  the  underside  of  the  plunger.  This  spring 
makes  a  sliding  electrical  contact  with  the  cavity  walls.  Not  shown  in  this  figure  are  two  rf  electrical 
feedthroughs  which  penetrate  about  1  nun  into  the  cavity  through  the  bottom  and  side  walls. 

Figure  4  (a)  shows  the  disassembled  rectangular  cavity  and  cylindrical  beam  tunnels.  The  tunnel 
diameter  is  0.33"  to  provide  good  cutoff  properties  at  12  GHz  and  to  reduce  space  charge  effects  in  the 
small  diameter  electron  beam.  In  order  to  make  the  cavity  comers  as  sharp  as  ]X>ssible,  the  cavity  was 
machined  in  four  pieces,  and  the  end  walls  were  brazed  to  the  cavity  body.  To  keep  the  beam  tunnels 
p>erpendicular  to  the  cavity  axis,  brazing  was  also  used  when  attaching  the  tunnels.  In  all  cases,  a 
nonmagnetic  Silcoro  75  braze  material  was  used.  In  the  assembled  cavity,  the  deviations  of  electron 
beam/drift  tunnel  wall  separation  are  less  than  1  mil.  The  assembled  cavity  and  cold  test  setup  is 
depicted  in  Rgs.  4(b)  and  4(c).  The  bronze  disk  on  the  cavity  cap  is  the  worm  gear,  with  the  tuning  shaft 
exiting  the  picture  on  either  side.  Microwave  power  is  fed  in  through  the  probe  in  the  cavity  side  wall. 
The  cavity  was  leak  checked  and  easily  obtained  pressures  of  10'^  torr. 

B.  Cold  Test  Results 

The  degree  of  coupling  of  coaxial  probes  into  this  rectangular  cavity  was  determined  via  the  High 
Frequency  Structure  Simulator  (HFSS)  in  collaboration  with  Code  6840  of  the  Naval  Research  Laboratory. 
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Figure  5  shows  the  CAD  geometry  with  the  rectangular  cavity  and  a  cylindrical  beam  tunnel  on  each  end. 
Coaxial  lines  couple  to  the  cavity  from  adjacent  sides.  In  the  figure  shown,  a  signal,  near  the  resonant 
frequency  of  11.80  GHz.  The  probe  loading  causes  the  decrease  of  the  resonant  frequency.  Tl.e  center 
conductor  of  the  coaxial  cables  penetrates  into  the  cavity  a  distance  .122  cm,  which  is  about  one-tenth  of 
the  way  across  the  cavity  and  is  flush  with  the  edge  of  the  beam  tunnel.  The  figure  shows  the  electric 
field  intensity  in  the  z-x  plane.  The  darker  the  color,  the  weaker  the  electric  field.  The  signal  is  supplied 
through  the  top  coaxial  cable.  Several  features  are  immediately  apparent.  First,  the  field  strength  in  the 
beam  tunnels  is  very  weak.  This  is  because  the  11.8  GHz  signal  is  so  far  below  cutoff.  Also  the  field 
intensity  in  the  coaxial  feed  is  quite  weak.  Thus  the  match  between  the  feed  and  the  cavity  is  good. 
Finally,  from  the  light  color  near  cavity  center,  we  see  a  strong  field  excited  in  the  appropriate  mode  in 
the  cavity.  Cross-coupling  between  the  coaxial  feeds  was  also  measured.  This  is  important  because  it  is 
we  want  the  probes  to  couple  to  only  one  mode.  (Experimentally  this  will  be  a  good  means  of  mode 
identification).  It  is  found  that  there  is  better  than  30  dB  isolation  between  the  probes  (signal  launched 
through  one  probe  and  picked  up  by  the  other).  When  a  signal  of  a  given  polarization  is  launched 
through  the  cutoff  tunnel,  similar  results  are  obtained. 

The  cavity  has  been  constructed  and  a  cold  test  of  cavity  resonant  frequency  and  (^-factor  has 
been  performed.  The  setup  is  pictured  in  Figs.  4(b)  and  4(c).  Both  reflection  and  transmission 
measurements  were  made  on  the  cavity.  The  resonant  frequencies  and  Q  values  were  found  for  the 
lowest  order  odd  modes  by  reflection  measurements.  The  frequencies  were  measured  using  a  sweep 
oscillator  and  cavity  wavemeter.  The  Q's  were  calculated  using  the  FWHM  of  the  resonance  curves. 
Only  odd  modes  can  be  detected  since  the  coupling  probe  is  at  the  cavity  center.  The  Q  values  and 
resonant  frequencies  are  given  in  the  table  below.  Note  that  the  resonant  frequency  of  the  TEjoi  mode. 


Resonant  Frequency 

Q 

^101 

11.810  GHz 

408 

^011 

12.062  GHz 

670 

^^013 

13.030  GHz 

434 

"^^015 

14.871  GHz 

708 

the  one  not  affected  by  the  tuning  plunger,  has  a  resonant  frequency  which  is  very  close  to  that  calculated 
by  HFSS.  The  quality  factors  are  quite  low,  also  in  agreement  with  HFSS.  By  adding  a  shorted  stub  in 
parallel  with  the  input  feeds,  these  Q's  can  be  increased  substantially.  A  transmission  measurement,  at 
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the  resonant  frequency  of  the  TEqh  mode,  indicates  that  there  is  about  24  dB  isolation  between  the  probes 
(again  in  reasonable  agreement  with  HFSS).  This  is  adequate  to  selectively  couple  the  two  modes. 
Finally,  the  resonances  were  unaffected  by  the  introduction  of  metal  at  the  drift  tube  ends.  This  indicates 
that  the  tubes  are  well  cutoff  to  the  radiation  frequencies  used  in  the  experiment. 


IV.  TUBE  DESIGN  AND  EXPERIMENTAL  INFRASTRUCTTURE 


A.  Tube  Design 

One  of  the  first  experimental  components  built  was  a  valve/bellows  assembly  to  be  attached  to 
the  electron  gun.  This  purpose  of  this  assembly  is  to  allow  the  electron  gun  to  stay  under  vacuum  while 
changes  in  the  nucrowave  circuit  are  made.  This  assembly  is  sketched  in  Fig.  6  and  a  photograph  is 
shown  in  Fig.  7.  The  gate  valve  is  straight-through  and  is  sealed  by  a  Viton  O-ring.  This  prevents  baking 
of  the  valve  at  high  temperature.  Good  vacuum  was  obtained  by  removing  the  Oring  and  baking  at  high 
temperature,  then  back-filling  the  system  with  dry  nitrogen  and  quickly  replacing  the  Oring.  The 
bellows  is  guided  by  two  control  rods  and  a  screw  rod  is  used  to  expand  the  bellows  against  the  force  of 
vacuum.  The  bellows  houses  a  beam  tube  which  can  be  inserted  into  the  valve  when  the  gate  is  open. 
This  prevents  the  electron  beam  from  being  exposed  to  the  non  uniform  internal  structure  of  the  gate 
valve. 

The  Varian  8063  gun  was  activated  by  applying  heater  power  and  a  small  dc  voltage  across  the 
mod-anode  /  cathode  gap.  The  power  supply  was  capable  of  supplying  10  nvA  of  current,  hence  we 
could  not  run  the  gun  to  veiy  high  voltage  (only  about  50  -1  00  V  dc).  The  scaling  looked  acceptable  and 
we  expect  full  current  at  high  voltage  operation. 

Capacitive  probes  are  used  as  beam  diagnostics  in  this  experiment.  The  capacitive  probe 
measures  directly  the  electron  beam  line  charge  density.  With  the  knowledge  of  total  beam  current,  the 
electron  axial  velocity  can  be  deduced.  The  important  parameter  a  (a  =  can  then  be  determined  by  a 

knowledge  of  the  beam  energy.  The  probe  assembly  consists  of  two  sets  of  SOD  SMA  vacuum 
feedthroughs  on  a  ConFlat  flange  as  depicted  in  Fig.  8.  The  feedthroughs  are  separated  by  180°  on  the 
circumference.  A  small  metal  disc  is  welded  on  the  tip  of  each  feedthrough  center  conductor  to  increase 
the  probe  capacitance.  The  discs  are  flush  with  the  drift  tube  inner  surface  and  electrically  isolated  from 
the  tube.  The  probe  signal  is  linearly  proportional  to  the  electron  beam  line  charge  density  when 
operated  in  a  high  impedance  mode  (1  MQ.)  Equation  (6)  shows  the  relationship  between  the  probe 
voltage  Vp  and  the  beam  line  charge  density, 

X  =  27lc(4-)Vp  (6) 

Ap 
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where  c  is  the  total  capacitance  of  the  measurement  set-up,  b'  is  the  effective  drift  tube  radius  and  Ap  is 
the  effective  probe  area.  The  constant  c(-^)  can  be  determined  experimentally  through  the  calibration 


process. 

In  the  calibration  set-up,  the  actual  electron  beam  was  replaced  by  a  metal  tube  with  a  similar 
radius  as  shown  in  Photograph  7(b).  The  metal  tube  was  pulsed  to  a  high  voltage,  V,  to  send  a 
displacement  current  to  the  probes.  The  relationship  between  V  and  Vp  can  be  written  as 


V=M^(i2!_)Vp  (7) 

^  Ap 

where  b  and  a  are  the  radii  for  the  drift  tube  and  the  center  metal  tube.  For  a  given  set-up  Equation  (7) 
uniquely  determines  the  constant  c(-^)  .  This  constant  then  can  be  used  to  determine  the  line  charge 

Ap 

density  through  Equation  (6).  The  axial  velocity  can  be  determined  as  V//  =  where  I  is  the  total  beam 


current. 

Fig.  9  shows  the  calibration  results  with  two  different  cable  length;  5  feet  and  13.1  feet.  One  of  the 
probes  is  found  to  be  slightly  misaligned  after  the  fabrication  and  is  believed  to  be  the  reason  for  the 
difference  in  probe  response.  Since  different  cable  lengths  represent  different  total  capacitance,  the  slopes 
are  therefore  different.  For  example  for  probe  #I  with  13.1  feet  cable,  the  constant  is  found  to  be 


7.14x10'®.  The  following  table  gives  the  values  for  the  constant  c(-^) 


with  various  cable  lengths. 


c(-*4-) 

Ap 

Probe  #1, 5  feet  cable 

3.01x10'® 

Probe  #1, 13.1  feet  cable 

7.14x10'® 

Probe  #2, 5  feet  cable 

3.84x10*® 

Probe  #2, 13.1  feet  cable 

9.09x10'® 

The  design  of  the  electron  beam  collector  was  closely  guided  by  the  POISSON  and  MAGIC  PIC 
simulation.  Fig.  10  shows  the  final  design  of  the  collector  including  the  magnet  polepiece.  In  order  to 
provide  the  flexibility  to  move  the  collector  in  the  axial  direction,  a  cup-shape  polepiece  is  incorporated. 
The  function  of  the  polepiece  is  to  create  an  abruptly  divergent  magnetic  field  at  the  entrance  of  the 
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collector  region  so  that  electrons  can  be  diverted  toward  the  collector  wall  where  cooling  is  provided. 
The  POISSON  simulation  generates  a  magnetic  field  distribution  according  to  the  electromagnet 
manufacturer's  specifications  and  the  proposed  polepiece  design.  This  simulation  result  is  then 
transferred  to  the  MAGIC  PIC  code  as  the  prescribed  field  for  electron  trajectory  simulation.  The 
resultant  electron  spatial  distribution  wilt  then  give  the  information  on  the  efficacy  of  the  overall  collector 
design.  The  criteria  for  success  is  that  after  passing  through  the  polepiece  the  electrons  will  deposit  their 
kinetic  energy  on  the  collector  body  where  a  proper  cooling  is  provided.  The  polepiece  thickness  and  the 
size  of  the  bore  hole  play  an  important  role  in  the  electrons  behavior  inside  the  collector.  The  iteration 
between  varying  the  pole  piece  shape  and  checking  the  electron  trajectories  is  continued  until  a  proper 
design  is  achieved.  Figure  11  shows  the  simulated  electron  trajectory  with  MAGIC  PIC  code  for  a 
polepiece  which  is  0.8"  thick  with  a  bore  hole  of  1.23"  in  diameter.  .  The  electrons  are  observed  to  be 
distributed  over  a  broad  area  on  the  collector  wall  (~  100  cm^.)  The  power  density  can  be  estimated  to  be 
«4.0kW/cm2. 

In  order  to  simultaneously  achieve  high  vacuum  capability  and  good  mechanical  stability,  a 
ceramic  ring  structure  was  chosen  to  provide  the  necessary  electrical  isolation  of  the  collector  from  the 
tube  body.  This  d.c.  break  involved  brazing  alununa  rings  to  30  mil  thick  Kovar  disks.  These  disks  were 
in  turn  electron  beam  welded  to  stainless  steel  flanges.  While  Kovar  was  needed  to  provide  a  good 
brazing  interface,  the  mass  of  the  material  was  minimized.  This  in  turn  minimized  the  possible 
perturbation  of  the  magnetic  field  in  the  collector  entrance.  The  insert  in  Fig.  10  depicts  the  actual 
dimenaon  of  the  structure. 

The  collector  body  was  made  of  304  stainless  steel  tube  with  3.0"  OD  and  0.1"  wall  thickness. 
Collector  cooling  is  provided  with  a  water  jacket  consisted  of  7-tum  1  /4"  copper  tubes.  With  a  cooling 
water  flow  rate  of  0.1  gpm,  the  collector  can  handle  a  beam  power  of  400  W  continuously. 

B.  Experimental  Infrastructure 

The  electromagnet  (manufactured  by  Premiere  Microwave,  Palo  Alto,  CA)  has  been  tested  and 
the  axial  magnetic  field  measured  with  an  F.  W.  Bell  9900  Gaussmeler  and  FAA  99(X)915  probe.  The 
electromagnet  was  powered  with  four  30  kW  DC  power  supplies.  In  the  measurement  shown  in  Fig.  13, 
DC  current  from  each  power  supply  was  80  A  which  generated  a  peak  magnetic  field  «  3.2  kG.  The 
current  was  measured  with  resistive  shunts.  The  magnetic  field  was  measured  at  a  05"  axial  increment 
throughout  the  whole  electromagnet  and  compared  with  POISSON  simulation.  Fig.  13  shows  both  the 
measurement  and  the  simulation  results.  It  was  observed  that  some  discrepancy  existed  in  region  at  15" 
away  from  one  end  of  the  magnet.  The  variation  was  about  2%  of  the  peak  magnitude.  The  actual  cause 
is  not  yet  clearly  understood.  The  electromagnet  has  also  been  running  with  90%  power  for  »  5  hours  at  a 
cooling  water  flow  rate  of  20  gpm  No  thermal  problems  were  observed. 
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The  main  purpose  of  the  trim  coil  is  that  in  conjunction  with  the  7  kG  solenoid,  it  produce  the 
desired  field  profile  near  the  electron  gun  cathode.  Specifically,  for  the  gun  that  is  being  used  in  this 
experiment,  it  has  been  determined  that  axial  magnetic  field  at  the  cathode  has  to  be  in  the  range  325-375 
Gauss,  and  that  the  field  be  relatively  constant  over  the  surface  of  the  cathode.  Other  issues  in  this  design 
are  that  the  B  field  has  to  increase  from  its  value  of  -350  Gauss  at  the  cathode  to  a  value  of  about  1  kG  in  a 
distance  of  1.5"  toward  the  solenoid.  This  requirement  comes  about  because  the  gun  beam  channel  is 
approximately  1.5"  away  from  the  cathode  center,  and  has  a  radius  of  about  0.05".  The  design  was 
achieved  by  using  the  POISSON  code.  Figure  14  shows  the  output  from  the  POISSON  simulation.  The 
trim  coil  is  run  at  20000  Ampere-turns,  and  its  right  edge  is  5"  away  from  the  left  edge  of  the  solenoid. 
The  solenoid  is  run  at  177  amps,  to  produce  a  axial  B  field  of  7  KG.  There  is  a  1.6"  pole  piece  at  the  edge  of 
the  solenoid,  and  the  trim  coil  is  surrounded  on  both  sides  by  iron  plates  of  0.25"  thickness,  and  with  an 
inner  radius  of  3.2"  and  outer  radius  of  9.5".  The  fields  are  me.  '  .ed  at  0.225"  off  the  axis,  since  the  center 
of  the  emitting  surface  of  the  cathode  is  also  0.225"  off  the  axis.  As  seen  from  Fig.  14,  a  relatively  flat  field 
of  about  350  Gauss  is  obtained  at  a  distance  of  about  2.5  inches  away  from  the  left  edge  of  the  solenoid 
with  this  set-up. 

The  trim  coil  (Applied  Magnetic  Products,  Oakland,  CA)  has  been  fabricated  and  delivered.  The 
coil  consists  of  5  pancake  sections  of  electromagnets.  A  copper  sheet  provides  heat  conduction  from  the 
center  of  each  section  toward  the  edge  where  it  is  then  was  removed  by  cooling  tubes.  The  trim  coil  is  5" 
long  with  a  center  bore  of  7".  It  is  powered  with  an  EMI  DC  power  supply.  Figure  15  shows  the 
comparison  between  measurement  and  POISSON  simulation.  The  result  indicates  that  the  measurement 
and  simulation  are  in  good  agreement.  A  trim  coil  support  which  is  capable  of  XYZ  motion  and 
adjustment  has  been  designed  and  manufactured.  In  principle  the  structure  is  capable  of  supporting  a 
weight  of  10(X)  lbs. 

Figure  16  shows  a  sketch  of  the  laboratory  setup.  The  lab  features  desk  space  with  two 
Macintosh  computers  and  an  IBM  486  (to  run  EGUN  and  MAGIC).  A  chemical  fume  hood  allows 
preparation  of  materials  for  a  high  vacuum  environment.  The  high  voltage  modulator,  equipment  racks, 
and  d.c.  p>ower  supplies  are  all  labeled.  The  magnets  rest  on  the  aluminum  optical  bench. 
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HGURE  CAPTIONS 


1:  Linear  coefficients  for  the  TEioi  mode  .  The  curves  represent  a'  (solid),  a"  (dash),  a'loi^ll 
(chaindot),  and  a"ioi,oil  (^’s)- 

2:  Phase  and  group  velocity  matching  and  separation  of  TEjoi  mode  from  the  next  higher  mode;  the 
TE102. 

3:  Machine  drawing  for  the  final  cavity  design. 

4:  Photographs  of  rectangular  cavity  and  cold  test  setup,  (a)  disassembled  cavity,  (b)  and  (c)  views  of 
cold  test. 

5:  HESS  analysis  of  electric  field  excitation  at  11.8  GHz  in  cavity.  Standing  wave  is  seen  in  the  feed  line 
and  evanscent  waves  are  present  in  beam  tunnels. 

6;  Sectional  drawing  of  the  electron  gun  and  valve/bellows  assembly. 

7:  Tube  components  (a)  electron  gun  and  valve/bellows  assembly,  (b)  calibration  setup  for  capacitive 
probes,  (c)  electron  beam  collector  with  water  cooling  coil. 

8:  Geometry  of  the  capacitive  probes  (a)  detail  on  vacuum  feedthrough  and  dimensions  of  probe  tip,  (b) 
mounting  configuration  for  probes  on  ConFlat  flange. 

9;  Capacitive  probe  response  for  cable  lengths  of  5  ft.  and  13.1  ft.  as  a  function  of  voltage  applied  to 
center  conductor. 

10:  Machine  drawing  for  electron  beam  collector  and  pole  piece. 

11:  MAGIC  simulation  of  10  keV  electron  beam  trajectories  into  the  collector  region. 

12:  General  layout  of  magnets  and  experiment  (a)  solenoidal  magnet  and  high  voltage  modulator,  (b) 
solenoidal  magnet  with  electron  gun  and  Lucite  gun  case  in  foreground  (c)  trim  coil  electromagnet 
and  px)wer  supply  with  solenoid  px)wer  supplies  in  background. 

13:  Magnetic  field  measurement  compjared  to  POISSON  simulation  with  80  A  in  coils  of  main  solenoid. 

14:  Superposition  of  trim  coil  field  and  solenoid  field  with  flat  region  at  cathode  pxjsition. 

15:  Magnetic  Peld  measurement  compared  to  POISSON  simulation  for  12  A  in  trim  coil. 

16:  Layout  of  laboratory. 
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